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SECTION  1 


INTRODUCTION 

l.i  BACKGROUND. 

Over  the  past  several  years,  the  Defense  Nuclear  Agency  has  funded 
the  Explosion  Effects  Division  (EED)  of  the  USAE  Waterways  Experiment 
Station  (WES)  to  develop  a  large-bore  (4-ft1  diameter)  vertical  gas  gun. 
The  primary  purpose  of  the  gas  gun  is  to  simulate  ground  shock  environ¬ 
ments  suitable  for  testing  ground  shock  measurement  devices.  A  poten¬ 
tial,  but  as  yet  unexplored,  secondary  purpose  for  the  gun  is  to  provide 
dynamic  material  properties  data  from  relatively  large  geologic  samples. 

Gas  guns  are  used  to  generate  localized  shock  environments  in 
materials,  for  shock  physics  studies,  by  high-velocity  impacts  of  a 
projectile  against  samples  of  the  materials  under  study.  A  gas  gun's 
ability  to  produce  controlled  and  repeatable  shock  inputs  is  an  attrac¬ 
tive  alternative  to  high  explosive  techniques,  which  are  more  commonly 
used  to  test  the  performance  of  ground  shock  transducers.  Typically,  gas 
guns  accelerate  projectiles  which  are  small  (6- in  dia.  or  less)  and 
achieve  high  velocities.  Small  gun  bores,  however,  limit  the  suitability 
of  existing  gas  guns  to  test  transducers  in  soils,  since  they  cannot  pro¬ 
duce  the  large  one -dimensional  fields  that  are  desirable  for  such  tests. 
The  purpose  of  the  large -bora  gun,  recently  constructed  by  WES,  is  to 
generate  these  large,  one -dimensional  stress  and  motion  fields  in  differ¬ 
ent  types  of  soils.  These  fields  can  then  be  used  for  controlled  tests 
of  ground  shock  transducers,  as  well  as  for  other  shock  physics  studies. 

Some  of  the  desired  attributes  of  the  4-ft  diameter  gun  are  that 
it: 

•  requires  no  explosive 

•  produces  repeatable  inputs 


‘A  table  of  factors  for  converting  non- SI  units  of  measurement  to  SI 
(metric)  units  in  presented  on  page  iv, 


•  provides  an  efficient  and  cost-effective  testing  r  athod 

•  generates  minimal  noise  levels  (allowing  use  in  relatively 
inhabited  areas) 

•  has  a  vertical  orientation  (for  easier  testing) 

1.2  DESCRIPTION  OF  4- FT  DIAMETER  VERTICAL  GAS  GUN. 

The  gun  was  designed  over  a  24 -month  period  beginning  in  1988 
(Ohrt,  (1988),  White,  (1990a),  White  (1990b)).  An  artist's  rendering  of 
the  4-ft  diameter  gas  gun  is  shown  in  Figure  1-1,  and  a  schematic  of  the 
gun  is  shown  in  Figure  1-2.  The  gun  consists  of  a  large  annular  pressure 
vessel  surrounding  a  vertical  barrel.  A  series  of  orifices  are  machined 
in  the  barrel  wall.  These  allow  the  compressed  air  from  the  vessel  to 
expand  into  the  barrel.  When  the  gun  io  in  the  cocked  position,  a 
projectile,  with  o-rings  placed  at  the  top  and  bottom,  straddles  the 
orifices,  and  prevents  the  compressed  air  from  being  released  into  the 
barrel.  The  projectile  is  held  in  place  >y  a  quick-release  trigger 
mechanism.  A  water  reaction  mass  fills  the  top  portion  of  the  barrel 
above  the  trigger  c’-^nism.  The  bottom  of  the  barrel  may  be  sealed  with 
a  diaphragm  to  allow  a  partial  vacuum  to  be  created  in  the  barrel  section 
below  the  projectile. 

To  fire  the  gun,  the  projectile  is  released.  The  weight  of  the 
projectile  causes  it  to  move  downward.  As  the  top  o-ring  clears  the 
orifices,  the  compressed  air  expands  into  the  barrel.  The  incoming  air 
simultaneously  drives  the  projectile  downward  and  the  water  reaction  mass 
upward.  Theoretically ,  the  mass  of  the  water  can  be  adjusted  so  that  the 
"bottom"  of  the  mass  of  water  will  exit  the  top  of  the  barrel  at  the  same 
time  that  the  projectile  clears  the  bottom  of  the  barrel. 

To  prove  that  this  working  concept  for  a  4-ft  diameter  vertical  gas 
gun  was  feasible,  a  12-in.  diameter  vertical  gas  gun  was  first  designed, 
built  (Joachim  (1985),  White  (1988)),  and  tested  (Ohrt  and  Welch  (1989)). 
The  results  of  the  dynamic  and  static  tests  conducted  with  the  12 -in.  gun 
indicated  that  the  4-ft  gun  concept  was  both  feasible  and  practical. 
Details  on  the  design  of  each  gun  and  the  testing  of  the  12-in.  gun,  are 
given  in  White,  et  al  (1991). 
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Flgurtt  \-l.  WES  A* ft  Ciijtetsv  V^xtlcal  Gas  Gun. 
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Figure  1-2.  Schotaatie  of  tl.o  4-ft  diiuseter  gas  gu:.. 
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1.3  SCOPE. 


This  report  presents  a  discussion  of  the  safety  tests  conducted  on 
the  gun.  Comparisons  between  predicted  and  measured  far  field  nuisance 
effects  (ground  motion  and  noise)  are  presented.  Also  presented  is  a 
mathematical  model  that  provides  upper  bound  estimates  for  the 
performance  of  the  4-ft  diameter  gas  gun,  i.e.,  the  projectile  velocity 
as  a  function  of  the  driving  vessel  pressure.  The  projectile  velocity 
measured  on  several  pressurized  tests  with  the  gun  are  compared  to  the 
predicted  values.  Conclusions  ars  drawn  on  the  results  of  the 
development  program  to  date. 


5 


SECTION  2 


SAFETY  TESTING 

2.1  HYDROSTATIC  TESTING. 

Hydrostatic  testing  (White,  1990c)  was  performed  to  ensure  the 
structural  integrity  of  the  pressure  vessel  and  the  projectile  o-ring 
seals  before  placing  the  4- ft  gas  gun  in  service.  The  initial  phase  of 
hydrostatic  testing  was  conducted  immediately  following  the  construction 
of  the  projectile.  The  second  phase  was  conducted  to  increase  the 
operating  level  to  the  maximum  design  level  of  2,070  kPa  (300  psi) . 

Figure  2-1  illustrates  the  test  set-up  for  both  the  initial 
hydrostatic  tests  and  vacuum  tests  (described  later).  The  gun  was 
oriented  horizontally  for  the  initial  tests.  At  the  top  of  the  pressure 
vessel,  two  high-pressure  female  couplings  were  attached,  one  as  a  fill 
line  and  the  other  as  a  bleed  line.  An  additional  coupling  was  attached 
at  the  bottom  of  the  tank  for  draining  the  water  after  the  hydrostatic 
tests  were  completed.  These  three  holes  are  used  during  normal  operation 
of  the  gas  gun  for  filling  the  tank  with  air,  attaching  a  relief  valve, 
and  attaching  a  pressure  gage. 

For  the  hydrostatic  tests,  the  vessel  was  first  filled  with  water. 
Compressed  air  was  then  applied  to  the  water-filled  tank.  The  only  air 
In  the  system  was  that  in  the  hose  connecting  the  compressor  to  the  gun. 
The  o-ring  material  used  for  sealing  around  the  projectile  for  the  first 
hydrostatic  test  was  a  70-durometer  nitrile  compound  with  a  19-aa 
(3/4-in.)  diameter.  Testing  revealed  that  this  o-ring  was  inadequate  for 
sealing  around  the  projectile.  The  failure  was  a  result  of  several 
lacerations  in  the  o-ring,  created  while  loading  the  projectile.  As  the 
o-ring  moved  along  the  orifices  in  the  barrel,  it  extruded  into  the 
orifices,  and  was  nicked  when  it  passed  beyond  the  top  of  the  orifices. 

The  second  hydrostatic  test  in  the  initial  series  used  16-aa 
(5/3-in.)  diameter  o- rings  (70  durooetcr  nitrile)  in  place  of  the  19-aa 
(3/4-in.)  projectile  o-rings.  The  smaller  thickness  o-ring  extruded  less 
into  the  orifices,  and  prevented  damage  to  the  o-ring.  The  vessel  was 
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Figure  2*1.  Test  setup  for  initial  hydrostatic  and  vacuum  testing 
of  the  4- ft  diameter  gas  gun. 

then  pressurized  to  2,275  kPa  (330  psi).  The  16-mm  o-ring  adequately 
sealed  around  the  projectile.  In  addition,  the  welds  of  the  pressure 
vessel  were  determined  to  be  sound,  verifying  the  structural  integrity  of 
the  vessel  for  conducting  tests  using  pressurized  air  up  to  1,520  kPa 
(220  psi).  A  third  hydrostatic  test,  identical  to  the  second  test, 
verified  these  results. 

A  second  series  of  hydrostatic  testa  was  performed  to  increase  the 
operating  level  of  the  gun  up  to  its  maximum  design  limit  of  2,070  kPa 
(300  psi).  These  tests  were  conducted  after  the  gun  had  been  placed  in 
use  and  tested  with  air  up  to  pressures  of  1,380  kPa  (200  psi).  The  set¬ 
up  for  the  second  series  of  hydrostatic  tests  is  shown  in  Figure  2-2. 

The  gun  was  in  its  normal  testing  orientation  (vertical)  for  these  tests. 
The  holes  in  the  vessel  that  are  normally  used  for  filling  the  tank  with 
air,  for  attachment  of  a  relief  valve,  and  for  attachment  of  a  pressure 
gage,  were  plugged  with  high-pressure  fittings.  Two  additional  holes 
were  cut  Into  the  vessel;  one  at  the  top  for  filling  the  vessel  with 
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water;  and  another  at  the  base  for  draining.  As  with  the  previous 
successful  hydrostatic  test,  16-sun  o-rings  were  used  to  seal  around  the 
projectile. 

After  filling  the  vessel  with  water,  nitrogen  was  used  to  increase 
pressure  within  the  pressure  vessel.  The  pressure  was  raised  to 
3,100  kP»i  (450  psi)  using  a  pressure  regulator  to  control  the  flow  of 
nitrogen  into  the  vessel.  A  minimal  loss  of  pressure  within  the  vessel 
was  noted  after  several  minutes  at  the  3,100  kPa  (450  psi)  level.  The 
test  was  repeated  to  verify  the  results. 

Undor  existing  Army  Corps  of  Engineers  safety  regulations,  tests 
can  be  conducted  with  the  gas  gun  using  pressurized  air  at  the  maximum 
design  limit  of  the  gun  (2,070  kPa,  300  psi);  i.e.,  two-thirds  of  the 
hydrostatic  test  pressure. 

2.2  VACUUM  TESTING. 

Vacuum  testing  (White,  1990c)  was  performed  by  evacuating  the 
barrel  of  the  gas  gun  beneath  the  projectile.  The  purpose  of  these  tests 
was  to  evaluate  (a)  the  performance  of  the  o-ring  seals  on  the  barrel, 

(b)  the  candidate  diaphragm  material  used  for  sealing  the  end  of  the 
barrel,  and  (c)  the  ability  of  the  trigger  mechanism  to  support  the 
vacuum  load  on  the  projectile. 

The  first  vacuum  test  used  19-ma  (3/4-in.;  diameter  o-rings,  made 
from  70-durometer  nitrile  material,  for  sealing  around  the  projectile.  A 
diaphragm  of  1-ma  (0. 040- in. )• thick  fiberglass  reinforced  polyester 
sealed  the  bottom  end  of  the  barrel.  After  running  the  vacuum  pump 
approximately  135  seconds,  the  diaphragm  failed  at  a  vacuum  level  of 
roughly  460  mm  of  mercury  (Hg). 

As  of  result  of  the  failure  of  the  1-mm  fiberglass  reinforced 
polyester,  a  second  test  was  conducted  using  a  12.7-ma  (1/2-in.)  aluminum 
plate  over  the  end  of  the  barrel.  This  test  used  projectile  o-rings  made 
from  16-mm  (5/8-in.)  diameter,  70-durometer  nitrile  material.  During 
this  test,  the  vacuum  pump  was  run  intermittently  while  observing  both 
the  vacuum  gage  and  the  deflection  of  the  diaphragm.  The  sealing 
capability  of  the  gun  was  assessed  by  observing  the  loss  of  vacuum  while 
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the  pump  was  off.  Approximately  25  mm  Hg  (3  percent  vacuum)  was  lost 
during  a  two-minute  span.  The  rate  of  loss  was  fairly  constant  over  the 
entire  range  of  vacuum. 

In  order  to  determine  the  time  required  to  pull  a  "full"  vacuum,  a 
third  test  was  conducted  during  which  the  pump  ran  continuously,  except 
for  brief  stops  for  taking  vacuum  level  readings.  The  maximum  attainable 
vacuum  was  approximately  740  mm  Hg,  which  was  achieved  in  15  minutes. 
After  running  the  pump  for  40  minutes,  no  change  in  the  vacuum  level  was 
evident.  The  pump  is  capable  of  pulling  down  to  a  vacuum  pressure  of 
10  microns  (759.958  mm  Hg),  but  there  are  sufficient  leaks  in  the  system 
that  prohibit  pulling  a  vacuum  greater  than  approximately  740  mm  Hg. 

The  rate  of  loss  of  vacuum  (12.5  mm  Hg/min)  was  determined  to  be 
acceptable.  It  is  not  a  requirement  to  hold  a  vacuum  for  an  extended 
period  of  time.  Because  of  the  displacement  capability  of  the  vacuum 
pump,  it  is  possible  to  simply  cycle  the  pump  on  and  off  Just  prior  to 
firing  the  gun  to  obtain  the  desired  vacuum  level,  up  to  740  mm  Hg. 

Also,  a  full  vacuum  is  not  required  to  fire  the  gun  (the  projectile  will 
begin  moving  within  the  barrel  under  its  own  weight,  once  released  by  the 
trigger  mechanism) . 

Also  observed  on  the  vacuum  test  was  movement  of  the  projectile. 

The  manner  in  which  the  quick-release  trigger  mechanism  is  constructed 
(Figure  2-3)  allows  for  a  slight  backward  rocking  motion  of  the  three 
supporting  latches  when  they  are  in  the  cocked  position  and  not  carrying 
a  load.  The  pretest  placement  of  the  projectile  was  such  that  the 
latches  of  the  trigger  mechanism  were  not  loaded.  Posttest,  the  latches 
were  tight  against  the  retaining  collar  and  could  not  be  rocked  backward, 
indicating  that  the  trigger  mechanism  was  Indeed  loaded  by  the  vacuum. 

In  future  tests,  the  only  loading  on  the  trigger  mechanism  (in  addition 
to  the  vacuum)  will  be  the  weight  of  the  projectile,  which  is  relatively 
small  when  compared  to  the  vacuum  loading;  i.e.,  14,500  N  (3,260  lbs)  as 
compared  to  155,600  N  (26,000  lbs). 

The  vacuum  tests  were  considered  successful.  The  critical 
components  of  the  trigger  mechanism  were  tested  and  found  to  perform 
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satisfactorily.  In  addition,  the  sealing  capability  of  the  o-rings 
within  the  gun  was  found  to  be  adequate. 

2.3  SAFETY  FEATURES  AND  SAFETY  PLAN. 

Safety  issues  have  received  the  highest  priority  in  the  development 
of  the  4-ft  diameter  gas  gun.  Some  of  the  safety  features  included  in 
the  design  and  siting  of  the  gas  gun  are: 

a.  Pressure  Vessel.  The  pressure  vessel  is  designed  for  a  maximum 
working  pressure  of  2,070  kPa  (300  psi).  This  design  includes  a  minimum 
factor  of  safety  of  four  in  all  components.  A  manway  is  incorporated  to 
provide  entry  to  the  vessel  for  periodic  inspection.  A  relief  valve  is 
attached  to  the  vessel  to  prevent  overpressurizing. 

b.  Trigger  Mechanism.  A  "potential  energy"  well  has  been 
incorporated  into  the  design  of  the  trigger  mechanism  (Figure  2-3)  used 
to  fire  the  gas  gun.  The  mating  surface  of  the  latch  and  the  retaining 
collar  of  the  trigger  mechanism  is  designed  such  that,  when  the  retaining 
collar  is  lifted  vertically,  the  latch  must  rotate  Into  its  housing. 

This  rotation  slightly  lifts  the  projectile  prior  to  firing.  Requiring 
the  application  of  a  significant  force  in  this  manner  in  order  to  fire 
the  gun  greatly  reduces  the  chance  of  a  misfire. 

c.  Operation  Controls .  The  controls  for  operation  of  the  gas  gun 
(Figure  2-4)  are  located  away  from  the  immediate  vicli.Iey  of  the  gun. 
These  controls  operate  the  "fill  system"  used  in  conducting  a  gas  gun 
tost.  The  fill  system  (illustrated  in  Figure  2-5)  is  comprised  of  three 
components;  the  pressure  system,  the  vacuum  system,  and  the  firing 
system.  The  pressure  and  vacuum  systems  incorporate  valves  to  bleed  off 
the  pressure  In  the  vessel  and  the  vacuum  in  the  barrel,  should  a  teat  be 
aborted,  or  in  the  event  of  a  misfire.  Also  included  in  the  control 
panel  is  a  master  switch  that  prevents  firing  the  gun  until  all  fill 
system  components  are  in  their  proper  state;  i.e.,  all  valves  closod  and 
the  air  compressor  and  vacuum  pump  are  shut  off.  The  Standard  Operating 
Procedure  (SOP)  for  conducting  a  test  is  posted  on  the  control  panel  in 
plain  view  of  the  operator. 
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Figure  2-5.  Fill  syatea  used  for  conducting  tests  with  the  gas  gun 


systematic  method  of  conducting  tests  in  a  safe  manner.  The  plan 
prescribes  the  safety  policies  and  procedures  for  testing  with  the  gas 
gun  and  applies  to  all  personnel  p.-ticipating  on  a  test.  The 
requirements  listed  below  were  mandatory  for  all  gas  gun  tests. 

a.  The  Standard  Operating  Procedure  is  followed  when  conducting  a 

test. 

b.  All  components  of  the  gun  are  inspected  on  a  regular  basis. 

c.  tohen  the  projectile  was  in  the  loaded  position  between  tests, 
the  battery  power  to  the  four -way  directional  valve  (used  for  firing  the 
gun)  is  disconnected.  Two  safety  chair*!  are  used  to  secure  the 
projectile's  position  within  the  barrel  between  tests. 

d.  No  personnel  are  allowed  in  close  proximity  (53  m  (175  ft))  to 
the  gun  when  a  test  is  bring  conducted. 

e.  A  warning  he  is  sounded  one  minute  prior  to  firing  the  gun. 

A  second  warning  is  sounded  thirty  seconds  prior  to  firing  the  gun.  Three 
short  blasts  are  joundad  to  count  down  the  final  seconds  prior  to  firing. 

f.  A  manned  road  block  is  placed  at  the  entrance  of  the  service 
road  to  the  test  s*te  prior  to  presurlzing  the  vessel. 

g.  A  standard  procedure  for  reentering  the  test  site  is  followed. 

As  a  minimum,  this  procedure  includes  a  visual  inspection  of  the  concrete 
pads,  the  platform  supposing  the  gun,  and  the  trench  shoring. 

h.  A  pump  is  maintained  on-site  to  l^move  any  water  in  the  trench. 
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SECTION  3 


NUISANCE  EFFECTS:  PREDICTIONS  AND  MEASUREMENTS 

There  was  no  experience  bare  available  for  estimating  the  probable 
range  of  nuisance  effects  associated  with  the  operation  of  a  gas  gun  as 
large  as  the  WES  4- ft  diameter  gun.  Accordingly,  an  asessment  had  to  be 
made  of  the  potential  nuisance -level  noise  and  far-field  ground  motions 
that  would  be  generated  by  the  operation  of  the  gun.  This  assessment 
included  the  development  of  upper  bound  predictive  equations  for  the 
induced  far-field  ground  motion,  and  careful  monitoring  of  the  induced 
noise  and  far-field  ground  motions,  during  the  test  to  insure  that  these 
quantities  were  within  acceptable  limits.  Predictive  equations  and  far- 
field  data  gathered  during  the  tests  are  presented  in  the  following 
sections . 


3.1  GROUND  MOTION. 


In  order  to  evaluate  potential  hazards  to  buildings  in  the  vicinity 
of  the  4-ft  gas  gun,  calculations  of  far-field  ground  motions  were 
performed  for  various  projectile  impact  velocities. 


Wallace  and  Fowler  (1973)  developed  a  relationship  for  peak 
vertical  particle  velocity  as  a  function  of  range  from  the  impacts  of 
spheres  dropped  onto  soil  surfaces  from  a  given  height.  The  equations 
convert  the  kinetic  energy  of  the  sphere's  impact  into  an  equivalent- 
explosion  energy  yield  for  TNT  explosives.  This  relationship  is 


Vt  «  1524 


(3.1) 


where , 

V,  -  peak  vertical  particle  velocity  for  impacts  (mo/see) 

E  **  yield  energy  (lbs  of  TNT) 

R  -  range  from  the  impact  or  explosion  point  (ft) 

By  determining  the  explosion  energy  yield,  in  pounds  of  TNT, 
equivalent  to  the  kinetic  energy  of  the  4-ft  gas  gun  projectile  at  the 
time  of  impact,  the  far-field  motions  can  be  predicted. 
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Figure  3-1  shows  the  predicted  upper  bound  impact  velocity  (as  a 
function  of  vessel  pressure)  for  a  1,478  kg  (3,260  lb)  projectile.  In 
this  calculation,  adiabatic  expansion  of  the  air  is  assumed,  and  the  mass 
of  the  air  and  friction  between  the  projectile  and  the  barrel  walls  is 
ignored.  A  detailed  development  of  the  mathematical  model  used  to 
predict  projectile  impact  velocity  is  presented  in  Section  4. 

The  kinetic  energy,  KE  (N-m),  of  the  projectile  at  the  time  of 
impact  is  given  by 

KE  *  YWj,V^  ^3‘2) 


where , 

Mp  -  projectile  mass  (kg) 

Vp  -  projectile  Impact  velocity  (m/sec) 

Substituting  Equation  (3.2)  into  Equation  (3.1),  and  using  a  TNT 
equivalence  of  1.912  x  104  N-m/lb-TNT  yields 


Vi  «  15240 


Mpvi 


N 


2K}(1.912  xlO*) 


(3.3) 


where , 

V,  -  peak  vertical  particle  velocity  (mm/soc) 

R  -  range  from  the  detonation  (ft) 

The  computer  code  FARTVEL  (listed  in  Appendix  C)  was  written  using 
Equation  (3.3)  to  calculate  the  peak  vertical  particle  velocity  for 
projectile  velocities  between  0  and  90  m/sec  (295  ft/sec)  for  various 
ranges  of  interest.  The  rangas  of  interest  are  listed  in  Table  3-1. 
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VESSEL  PRESSURE  (KPa) 


Figure  3-1.  Predicted  velocity  for  a  1,478-kg  (3,260-lb)  projectile 
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Table  3-1.  Ranges  of  interest  for  testing  with 
the  4-ft  gun. 

Location 

Range,  m  (ft) 

Nearest  (uninhabited)  bldg 

9  (30) 

Control  Trailer,  location 
of  nearest  personnel 

53  (175) 

Nearest  inhabited  WES  bldg, 
not  associated  with  testing 

78  (255) 

Guard  at  entrance  road  to 
site 

100  (325) 

Nearest  residential  bldg 
(outside  station  boundary) 

233  (765) 

The  threshold  of  human  perception  of  ground  vibration  is 
significantly  lower  than  the  levels  associated  with  the  onset  of 
structural  damage.  Subjective  human  response  to  vibratory  ground  motion, 
based  on  earthquake  studies,  has  shown  that  motions  of  0.1  mm/sec 
amplitude  are  the  absolute  lower  limit  of  human  perception,  and 
amplitudes  of  less  than  1  mm/sec  are  rarely  perceived  for  short-period, 
explosion-produced  motions  (Ristvet,  1987).  Ristvet  lists  a  level  of 
20.3  mm/sec  as  "unpleasant, "  and  Siskind,  et  al  (1990)  gives  17.8  mm/sec 
as  the  level  of  "discomfort,"  or  producing  a  "startle"  effect.  Siskind 
also  lists  thresholds  of  56  mm/sec  and  112  mm/sec  for  an  onset  of 
interference  with  activity  or  proficiency,  and  a  health  limit, 
respectively. 

Peak  particle  velocity  is  usually  taken  as  the  significant 
parameter  in  the  development  of  damage  criteria  for  structures.  Listed 
in  Table  3-2  is  a  summary  of  damage  thresholds  for  residential 
structures,  taken  from  several  references.  In  general,  these  criteria 
state  that  no  structural  damage  should  occur  below  a  peak  particle 
velocity  of  50  mm/sec. 

The  results  of  PARTVEL  calculations  for  the  9-m  (30-ft)  range  are 
shown  in  Figure  3-2.  Though  this  prediction  indicates  that  a  building  at 
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Table  3-2. 

Structural  damage  thresholds  from  references. 

Particle  Velocity  Damage  Threshold 

Damage  Type 

Langefors,  et  al 
(1958) 

Nicholls,  et  al 
(1971) 

McPherson 

(1989) 

None 

<  71  mm/sec 
(2.8  in. /sec) 

<  50  mm/sec 
(2.0  in. /sec) 

<  50  mm/sec 
(2.0  in. /sec) 

Fine  plaster 
cracks 

109  mm/sec 
(4.3  in. /sec) 

50-100  mm/sec 
(2-4  in. /sec) 

- 

Plaster  and 
masonry  wall 
cracking/  minor 
structure 

160  mm/sec 
(6.3  in. /sec) 

100-178  mm/sec 
(4-7  in. /sec) 

137  mm/sec 
(5.4  in. /sec) 

Major 
structural 
damage/  serious 
cracking 

231  mm/sec 
(9.1  in. /sec) 

>  178  mm/sec 
(7  in. /sec) 

193  mm/sec 
(7.6  in. /sec) 

this  range  could  incur  structural  damage  when  testing  with  projectile 
velocities  greater  than  27  m/sec  (90  ft/sec) ,  no  damage  to  the  storage 
building,  located  at  the  9-m  range,  was  noticed  after  testing  at  even  the 
highest  projectile  velocities. 

Shown  in  Figure  3-3  are  the  predicted  peak  vertical  particle 
velocities  as  a  function  of  projectile  velocity  for  the  remaining  ranges 
of  interest.  Triaxial  motion  measurements  (vertical,  radial  and 
tangential)  were  made  at  three  far-field  ranges  on  several  tests  with  the 
gas  gun.  These  ranges  were  53  m,  78  m,  and  233  m.  The  peak  value  of  the 
three  measurements,  at  a  given  range,  for  a  single  gas  gun  test  is  shown 
in  Figure  3-3.  This  data  is  tabulated  in  Table  3-3.  Note  that,  at  the 
nearest  residential  building  outside  the  WES  boundary  (233  m) ,  the 
measured  values  of  ground  motion  are  two  orders  of  magnitude  lass  than 
that  level  typically  regarded  es  &  hazard  to  residential  structures.  At 
the  location  of  the  nearest  personnel  (the  control  trailer  at  53  a),  the 
levels  are  a  factor  of  five  lower  than  the  allowable  threshold 
(50  mm/sec,  or  2  in. /sec).  Personnel  located  at  the  53-,  78-,  and  100-m 
ranges  could  feel  the  impact  for  most  projectile  velocities. 
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PROJECTILE  VELOCITY  (M/SEC) 


Figure  3*2,  Predicted  peek  particle  velocity  as  a  function  of 
projectile  velocity  for  the  9«a  (30-ft)  range. 
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PROJECTILE  VELOCITY  (M/SEC) 


Figure  3-3.  Predicted  peak  particle  velocity  as  a  function  of 
projectile  velocity  for  various  far-field  ranges. 
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Max  velocity  from  vertical  (v) ,  radial  (r)  ,  or  tangential  (t)  component 
Velocity  based  on  one  set  of  T0A  pins 
88%  of  predicted  value 
)  —  Ho  Data 
?  -  Not  Fielded 


Table  3-3.  Sound  pressure  level,  far- field  velocity,  and  airblast  data  from  WES 
4- ft  Disaster  Vertical  Gas  Gun  tests.  (Continued) 


The  calculations  presented  in  Figures  3-2  and  3-3  do  not  consider 
the  topography  of  the  area,  which  perhaps  accounts  for  the  discrepancy 
between  predicted  and  measured  values.  The  elevation  of  the  gun  is 
50.3  m  (165.1  ft).  An  embankment  adjacent  to  the  gun  rises  to  an 
elevation  of  54.9  m  (180  ft)  before  dropping  to  41.6  m  (136.6  ft)  at 
Durden  Creek,  some  143  m  (470  ft)  away.  The  elevation  rises  to  57.7  m 
(189.2  ft)  at  the  nearest  residential  structure  located  outside  of  the 
WES  boundary.  These  severe  changes  in  grade  that  exist  between  the  gun 
and  locations  of  interest  apparently  aided  in  attenuating  the  ground 
shock. 

3.2  AIRBLAST. 

During  a  gas  gun  test,  the  pressurized  air  driving  the  projectile 
and  reaction  mass  will  vent  at  both  the  top  and  bottom  of  the  barrel. 

This  pressure  release  into  the  atmosphere  will  cause  an  airblast  wave  in 
the  vicinity  of  the  gas  gun.  Personnel  in  the  area  will  sense  the 
airblast  wave  by  hearing  it,  and  if  close  enough,  by  feeling  it.  In 
order  to  evaluate  potential  hazards  to  personnel  (hearing  damage)  and 
buildings  (window  breakage),  calculations  were  performed  to  predict  the 
airblast/nuisance  levels  from  testing  with  the  gun.  Sound  pressure  level 
measurements  were  made  during  several  tests  at  all  five  ranges  of 
interest  (see  Table  3-1).  In  addition,  measurements  of  airblast  were 
made  at  three  far- field  ranges  (53,  78,  and  233  m)  on  several  tests. 

Sound  pressure  level  data  gathered  during  tests  with  the  WES  12 -in 
diameter  vertical  gas  gun  was  used  to  predict  the  airblast  level,  for  the 
ranges  of  interest,  when  testing  with  the  4-ft  gun.  The  noise  level 
associated  with  the  airblast  at  these  ranges  was  then  determined. 

During  the  performance  evaluation  tests  of  the  12-in  gun,  sound 
pressure  level  measurements  were  made  at  the  9-  and  33. 5 -m  ranges  for 
several  different  vessel  pressures.  This  data  is  presented  in 
Figure  3-4.  Note  from  this  figure  that,  with  the  exception  of  the  two 
measurements  made  at  the  33.5-m  (110-ft)  range  for  a  vessel  pressure  of 
345  kPa  (50  psi),  the  data  i^  relatively  constant  for  a  given  range. 
Therefore,  an  average  value  of  the  sound  pressure  level  (using  all  data 
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VESSEL  PRESSURE  (KPa) 


Figure  3 -4.  Sound  pressure  level  measurements  made  when  testing 
with  the  U&S  12-in.  Diameter  Vertical  Gas  Gun. 
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points)  was  determined  for  each  range.  The  peak  airblast  pressure 
amplitude  associated  with  a  given  sound  pressure  level  measurement  may  be 
determined  from  the  relations  (Kinsler,  et  al,  1982) 

p-vTp.  <3-4> 


and 


SPL  -  20 


(3.5) 


where , 

P  -  peak  pressure  amplitude 
P,  -  effective  pressure 
P*,  «■  reference  pressure  (20  iPa) 
SPL  -  sound  pressure  level  (dB) 


Solving  for  P4  in  Equation  (3.5)  and  substituting  into  Equation  (3.4) 
yields  the  peak  pressure  amplitude,  P  (in  Pascals): 

P  -  (3,6) 

Substituting  the  average  value  of  the  SPL  data  for  the  12- in  gun  into 
Equation  (3.6)  yields  the  peak  pressure  amplitude  at  the  two  ranges. 

In  order  to  use  the  12 -in.  gun  data  for  predicting  the  peak 
pressure  amplitude  at  various  ranges  of  interest  for  the  4-ft  gun,  each 
range  was  normalized  by  dividing  by  the  gun  diameter  (11.5  in.).  The 
results  of  Equation  (3.6)  for  the  12 -in.  gun  data  are  tabulated  as  a 
function  of  the  normalized  range  in  Table  3-4  and  presented  graphically 
in  Figure  3-5.  Listed  in  the  Table  3-5  are  the  peak  pressure  amplitude 
and  SPL  predictions  (scaled  from  the  linear  fit  to  the  12- in.  gun  data 
presented  In  Figure  3-5)  at  the  ranges  of  interest  for  the  4-ft  gun.  The 
ranges  are  also  expressed  in  normalized  values  (range/gun  diameter). 
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Table  3- 

4.  Airblast/nuisance  data  for  the  12 

-in  gun. 

Range 

Normalized  Range 

Peak  Pressure 

Avg.  SPL 

m  (ft) 

Pa  (psi) 

dB 

9  (30) 

31 

348  (0.05) 

141.8 

33.o  (110) 

115 

137  (0.02) 

133.7 

Table  3-5.  Predicted  airblast/nuisance  data  for  the  4-ft  gun. 


Range 

Normalized  Range 

Peak  Pressure 

SPL 

m  (ft) 

Pa  (psi) 

dB 

9  (30) 

7.5 

931  (0.135) 

150.3 

53  (175) 

43.75 

269  (0.039) 

139.6 

78  (255) 

63.75 

207  (0.030) 

137.3 

100  (325) 

81.25 

176  (0.0255) 

135.9 

233  (765) 

171.25 

105  (0.0152) 

131.4 

The  sound  pressure  level  corresponding  to  the  scaled  values  of 
pressure  for  the  4-ft  gun  way  be  determined  by  solving  for  P.  in 
Equation  (3.4)  and  substituting  into  Equation  (3.5). 


SPL  ■  ZOLod - £ - )  (3-7> 

\  (20S-6)  ) 

where  P  has  units  of  Pa  and  SPL  has  units  of  dfl.  By  substituting  the 
scaled  peak  pressure  amplitudes  for  the  4-ft  gun  into  Equation  (3.7),  the 
predicted  SPL's  can  be  tabulated  (Table  3-5)  and  presented  graphically 
(Figure  3-6). 

OSHA  regulations  state  that  exposure  to  impulsive  or  impact  noise 
Should  not  exceed  140  dB  peak  sound  pressure  level.  This  threshold  is 
included  in  Figure  3-6  for  comparison  to  predicted  and  measured  noise 
level®  at  the  various  ranges  of  interest  for  testing  with  the  4-ft  gun. 
For  the  4-ft  gun,  the  9-m  range  is  the  only  one  of  concern  that  lies 
above  the  140  dB  limit.  As  mentioned  previously,  the  closesC  range  at 
which  personnel  were  located  during  testing  was  53  o.  As  with  the  ground 
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PRESSURE  (Pa) 


Figure  3-5.  Airblast  pressure  measurements  and  predictions  as  a 
function  of  normalised  range  for  gas  gun  testing. 
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NORMALIZED  RANGE  (RANGE/GUN  DIAMETER) 


Figure  3-6.  Measured  and  predicted  sound  pressure  levels  as  a 
function  of  normalised  range  for  the  4-ft  gas  gun 
tests. 
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shock  calculations,  this  analysis  does  not  take  into  account  the  natural 
terrain  (hills  and  trees)  effects  that  will  tend  to  mitigate  the  airblast 
wave  and  thus  reduce  noise  levels. 

Presented  in  Figure  3-7  are  the  sound  pressure  level  data,  for  the 
various  ranges  of  interest,  as  a  function  of  vessel  pressure.  Note  from 
this  figure  and  Figure  3-6  that,  although  the  measured  levels  are 
relatively  high,  they  are  well  below  limits  requiring  safety  precautions 
(i.e.,  ear  protection).  Personnel  located  approximately  3.2  km  (2  miles) 
from  the  gun,  aware  that  a  test  was  imminent,  heard  the  gun  during  tests 
at  the  very  highest  levels;  i.e,,  1,725  kPa  (250  psi)  and  greater. 

The  most  likely  component  of  an  ordinary  structure  to  sustain 
damage  from  a  blast  wave  is  a  window.  Therefore  most  damage  criteria  are 
based  on  window  pane  breakage.  Not  all  window  panes  will  break  at  the 
same  blast  pressure  level,  so  the  study  of  this  phenomena  is 
statistically  based.  Some  large  plate  glass  windows  may  break  at  a 
pressure  level  of  207  Pa  (0,03  psi).  At  690  Pa  (0.1  psi)  some  windows 
break  and  at  6,900  Pa  (1.0  psi)  most  windows  break.  At  a  pressure  level 
of  20.6  kPa  (3.0  psi)  conventional  structures  are  severely  damaged 
(Blaster's  Handbook,  1980).  The  207  Pa  threshold  for  window  pane  damage 
is  included  in  Figure  3-5.  From  this  figure,  or  from  Table  3-5,  it  was 
determined  that  windows  in  buildings  located  at  the  9,  53,  and  78-m 
ranges  might  sustain  damage.  The  probability  of  window  damage  cannot  be 
determined  from  this  analysis.  However,  it  is  interesting  to  note  that 
the  close-in  data  point  for  the  12-in.  gun  is  located  above  the  207  Pa 
(0.03  psi)  threshold  and  windows  at  that  range  were  not  damaged  during 
performance  testing  with  the  12-in.  gun. 

The  windows  of  the  building  at  the  9-o  (30-ft)  range  were  covered 
with  a  sheet  of  plywood  to  provent  possible  breakage  during  a  tost  with 
the  A-£t  gun.  No  windows  were  damaged,  at  any  location,  during  any  test 
with  the  A*ft  gas  gun. 
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Figure  3-7.  Measured  sound  pressure  levels  generated  by 
the  4 -ft  gas  gun  as  a  function  of  vessel 
pressure,  for  various  ranges  of  interest. 
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SECTION  4 


MATHEMATICAL  MODEL 

A  mathematical  model  of  the  gas  gun  operation  was  formulated  by 
Welch  and  Ohrt  (White,  et  al  (1991))  in  order  to  calculate  projectile 
velocity  as  a  function  of  initial  vessel  pressure. 

Consider  the  simplified  view  of  the  gas  gun  shown  in  Figure  4*1. 

The  dashed  lines  show  the  initial  positions  of  the  projectile  and 
reaction  mass  (masses  M,  and  M,,  respectively).  The  trigger  mechanism  and 
its  supporting  members  are  located  in  the  chamber  between  the  two  masses. 
This  initial  volume  is  denoted  by  V^,.  and  the  volume  of  the  pressure 
vessel  by  V^.  The  area  of  the  barrel  is  given  by  A.  The  displacement  of 
the  projectile  is  given  by  the  variable  x  and  that  of  the  reaction  mass 
is  given  by  the  variable  y.  Hence,  the  total  volume  at  any  time  is 

v  -  ♦  v ^  *  AU>y)  (*•!> 

Let  the  pressure  at  any  time  after  firing  be  P.  Since  the  areas  of 
the  projectile  and  reaction  mass  are  the  same,  the  gas  exerts  an  equal 
force,  F  -  FA,  on  each.  Ignoring  frictional  forces  between  the  barrel 
walls  and  the  projectile  and  reaction  masses,  Newton's  second  law 
produces 

F  ■  Afpi  ■  Mrp 

Rearranging  and  integrating  twice  yields 

y  m  (4*2) 

W|> 

Substituting  Equation  (4.2)  into  Equation  (4.1)  yields 

.^1.  A)  (4.3) 

For  a  given  projectile  displacement,  x,  Equation  (4.3)  provide#  the 
volume  of  gas,  V,  of  the  system.  Assuming  adiabatic  expansion  of  the 
gas,  the  pressure,  P,  of  that  volume  has  the  relation 
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BEFORE  FIRING 


AFTER  FIRING 


t; 


REACTION  MASS 
-MASS,  Mr 
-AREA,  A 


PROJECTILE 
MASS,  Mp 
-AREA,  A 


Figure  4-1.  Parameters  used  in  developing  the  mathematical  model. 
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gas,  the  pressure,  ?,  of  that  volume  has  the  relation 


PV*  ■  constant 


where  y  is  the  ratio  of  specific  heats  for  the  gas  (about  1.4  for  air). 
If  the  pressure  and  volume  are  known  at  any  time,  say  Px  end  Vx,  then  the 
pressure  and  volume  as  a  function  of  time  is  given  by, 


Pit) 


PjVi'4 

(V(t))1-* 


(4.4) 


Before  the  projectile  and  reaction  masses  are  driven  by  the  gas, 
the  gas  must  first  expand  into  the  chamber  between  the  two  .masses.  The 
pressure  (Ppy)  and  volume  (Vpv)  of  the  pressure  vessel,  and  the  chamber 
volume  (VolMa)  are  known.  Using  Equation  (4.4),  the  initial  pressure 
acting  on  the  masses  can  be  calculated  as 


fiait 


P  V1,4 
rj>vvpy 


<V 


rchsmf 


1 1*4 


(4.5) 


The  pressure  acting  on  the  projectile  that  produces  a  given  projectile 
displacement,  x,  may  be  found  by  substituting  Equations  (4.3)  and  (4.5) 
into  Equation  (4.4), 


P  •  P 


ini  d 


V 


v-  4  •  -j) 


(4.6) 


The  work  done  on  the  projectile  by  this  pressure  is  given  by 

W  -  Jdfc'  -  )PctV  -  lo'PAdx 

where  Lp  is  the  distance  traveled  by  the  projectile.  Substituting 
Equation  (4.6)  in  the  last  of  these  expressions  yields 


fhp< , 

i" 

Jo  ifli< 

Vp,  *  ♦ 

Ik) 

**) 

1.4 

Adx 


(4.7) 


After  integrating,  this  expression  becomes 
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(4.8) 


2-5J»iait(Vov  ♦  W 

J 

1  - 

( _ 1 

0.4 

(-4) 

J. 

VPv  *  Vpfcu.  ♦ 

Mr)  J 

The  velocity  of  the  projectile,  Velp,  may  be  found  by  equating  the 
work  done  on  the  projectile  with  the  kinetic  energy  of  the  projectile. 
This  ignores  the  kinetic  energy  of  the  compressed  air  moving  behind  the 
projectile  and  is  an  upper  bound  calculation  of  the  projectile's 
velocity.  Thus, 

KB  »  -itfpVeip  ■  W 


Substituting  Equation  (4.8)  and  solving  for  Velp  yields, 


Vel, 


1  - 

0.4 

4-4) 

(4.9) 


Equation  (4.9)  gives  the  velocity  of  the  projectile,  for  a  given 
vessel  pressure,  after  it  has  travelled  a  distance  Lp.  Recall  that  P^t, 
in  this  equation  is  given  by  Equation  (4.5).  The  reaction  mass,  Mr,  is 
selected  using  Equation  (4.2)  by  setting  the  height  of  the  reaction  mass 
equal  to  the  distance  it  must  travel,  l*.  The  distance  the  projectile 
travels  before  leaving  the  barrel  is  Lpj  therefore,  the  travel  distance 
of  the  reaction  mass,  1*,  is  given  by 


4, 


"fit 

Mr  P 


(4.10) 


But  «  ALjP,  where  p  is  the  density  of  the  reaction  mass.  Substituting 
this  into  Equation  (4.10)  and  rearranging,  we  have, 


Ml 

Ap 
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Hence, 


Mz  «  <4*U> 

The  predicted  projectile  velocity  for  a  1,478  kg  (3,260  lb) 
projectile,  shown  in  Figure  3-1,  was  calculated  using  the  computer  code 
GG4PV,  listed  at  Appendix  D.  This  code  uses  Equations  (4.5)  and  (4.11) 
in  Equation  (4.9)  to  calculate  the  projectile  velocity  as  a  function  of 
vessel  pressure.  The  volume  of  the  pressure  vessel,  Vr,  is  constant 
(10.05  m5,  355  ft3)  as  is  the  area,  A,  of  the  barrel  (1.17  a*,  12.57  ft3). 
The  mass  of  the  projectile,  M,,  the  distance  traveled  by  the  projectile, 
ly,  and  the  chamber  volume,  V^,  vary  depending  of  the  configuration  of 
the  projectile,  The  value  for  each  is  typically  1478  kg  (3260  lb,), 

3.26  a  (10.7  ft),  1.8  n?  (63.7  ft3),  respectively. 
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SECTION  5 

PERFORMANCE  TESTS  OF  4-FT  GAS  GUN 
5.1  PROJECTILE  CONFIGURATION. 

The  projectile  design  currently  used  in  the  4-ft  gas  gun  is  shown 
in  Figure  5-1.  This  design  incorporates  the  initial  projectile  design 
(White,  1990b  and  White,  et  al,  1991)  by  using  it  as  a  carriage  for  an 
impact  plate.  Improvements  to  the  initial  design  of  the  projectile 
include : 

•  The  thickness  of  the  bottom  plate  of  the  carriage  was  increased 

•  The  thickness  of  the  material  used  for  the  stiffening  gussets 
inside  the  carriage  was  increased 

•  The  weight-reducing  holes  in  the  gussets  of  the  former  design 
were  eliminated 

•  Additional  smaller  gussets  were  added  inside  the  carriage 

«  Closed- cell  polyurethane  foam  was  added  between  the  carriage  and 
the  impact  plate. 

The  closed-cell  polyurethane  foam,  sandwiched  between  the  bottom 
plate  of  the  carriage  and  the  impact  plate,  is  used  to  limit  the  stress 
acting  on  the  bottom  of  the  carriage.  As  testing  levels  increase  (higher 
vessel  pressures  and  projectile  velocities),  thicker  pieces  of  foam  are 
required  to  absorb  the  impact  energy,  and  thus  protect  the  carriage  from 
damage.  Foam  thicknesses  varied  from  10  cm  (4  in.)  to  30  cm  (12  in.) 
during  recent  tests  with  the  gun.  The  mass  of  the  projectile  varied  very 
slightly  as  a  result  of  changing  the  foam  thickness.  All  tests  used  the 
projectile  configuration  of  Figure  5*1.  The  carriage  was  not  damaged 
during  the  test  series.  Two  5-cm  (2-in.)  thick  impact  plates  were  used 
during  the  test  series.  Other  than  occasionally  stripping  out  the 
threaded  holes  in  the  impact  plate  during  a  test,  the  plates  were  not 
damaged. 
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Figure  5-1.  Cross-section  of  the  gas  gun  projectile. 


5.2  PROJECTILE  VELOCITY  AND  PLANARITY  MEASUREMENTS. 


Six  piezoelectric  pins,  located  at  the  base  of  the  barrel  of  the 
gun,  were  used  to  monitor  the  velocity  and  planarity  of  the  projectile  as 
it  exited  the  barrel.  Shown  in  Figure  5-2  are  a  plan  view  and  a  section 
view  of  the  location  of  the  6  pin3.  Pins  T-1A  and  T-3A  are  located 
25.4  mm  (1  in.)  above  the  plane  containing  the  4  pins  designated  T-1B, 
T-2,  T-3B,  and  T-4.  A  schematic  showing  the  orientation  of  the 
piezoelectric  pin  within  its  bushing,  and  the  projectile  within  the 
barrel  of  the  gun,  is  shown  in  Figure  5-3.  A  Time -of -Arrival  Data  System 
(TOADS)  box,  running  at  a  clock  speed  of  1  MHz,  was  used  to  record  the 
time-of-arrival  (TOA)  of  the  projectile  at  each  pin  location.  Two 
measurements  of  projectile  velocity  were  computed  by  comparing  the  TOA  of 
the  projectile  at  the  locations  T-1A/T-1B  and  T-3A/T-3B.  The  planarity 
of  travel  of  the  projectile  was  computed  by  comparing  the  measurements  at 
pin  locations  T-1B,  T-2,  T-3B,  and  T-4. 

Several  factors  determine  the  accuracy  of  the  velocity /planarity 
measurements.  These  include  the  accuracy  of  machining  the  holes  in  the 
barrel  and  the  holes  in  the  pin  bushing,  the  flatness  of  the  impact 
plate,  the  assembly  of  the  projectile,  and  the  positioning  of  the  pins 
within  the  bushing.  Considering  these  factors,  the  accuracy  of  the 
location  of  each  pin  was  determined  to  be  ±  0.6  mm  (0.023  in.). 

Therefore,  the  actual  distance  between  two  pins  used  to  measure  the 
projeceilo  velocity  could  be  as  much  as  26.6  mm  (1.046  in.)  or  as  little 
as  24.2  mm  (0.954  in.).  Hence,  the  accuracy  of  the  velocity  measurement 
is  ±  4.6  percent.  Likewise,  pins  assumed  to  be  located  in  a  common 
horizontal  plane  could  have  a  relative  difference  in  their  vertical 
position  of  1.2  mm.  The  tolerance  on  the  planarity  measurement  depends 
on  the  location  of  the  two  pins  exhibiting  the  greatest  disparity  in  TOA. 
For  pins  located  90*  apart  the  tolerance  is  ±1.439  miiliradians  and  for 
pins  located  180’  apart  the  tolerance  is  ±1.018  miiliradians.  The 
frequency  response  of  the  piezoelectric  pin  and  the  recording  system  is 
high  enough  to  eliminate  any  significant  contribution  to  error  of  those 
measurements , 
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Figure  5-2.  Piezoelectric  pin  locations  for  velocity  and 
planar icy  aaasutenencs. 
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Figure  5*3.  Schematic  of  orientation  of  peisoelectric  pin  used 
in  measuring  projectile  velocity  and  planarity. 
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The  FORTRAN  program  TOACHECK,  listed  in  Appendix  E,  was  used  to 
evaluate  the  TOA  data  from  each  gas  gun  test.  This  program  calculates 
the  velocity  and  the  associated  accuracy  for  each  of  the  two  sets  of  pins 
T-1A/T-1B  and  T-3A/T-3B.  It  also  calculates  the  maximum  angle  of 
rotation  and  the  associated  accuracy  between  any  three  cf  the  four  pin 
locations  T-1B,  T-2,  T-3B,  and  T-4.  The  redundant  planarity  measurement 
is  used  to  assess  the  consistency  of  the  TOA  data.  A  plane  is  fit 
through  three  of  the  data  points  and  the  location  of  the  fourth  point  is 
projected.  The  measured  value  is  then  compared  to  this  calculated  value. 
This  procedure  is  conducted  for  each  possible  combination  (four  total). 
The  analysis  indicated  that  the  method  used  for  measuring  the  velocity 
and  planarity  of  the  projectile  is  consistent  and  congruous.  (Note:  a 
detailed  discussion  of  this  analysis  is  beyond  the  scope  of  this  report) . 

5.3  TEST  RESULTS. 

A  total  of  22  tests  were  conducted  with  the  WES  4- ft  Diameter 
Vertical  Gas  Gun  (see  Table  5-1  -  the  first  two  tests  were  preliminary 
exercises,  and  have  been  omitted  from  the  table).  The  vessel  pressure  in 
these  tests  varied  from  0  to  2,070  k?a  (300  psi)  and  the  projectile 
velocities  varied  from  5.4  to  7'i  m/sec  (17.6  to  232.8  ft/sec).  Multiple 
tests  were  conducted  at  each  of  several  vessel  pressure  levels  to 
investigate  the  repeatability  of  the  gun's  performance  with  regard  to 
nuisance  factors  (ground  vibration  and  noise)  associated  with  testing  and 
projectile  velocity. 

The  measured  projectile  velocities  for  tests  with  the  current  gas 
gun  projectile  are  shown  in  Figure  5-4.  The  circles  in  this  figure 
represent  predicted  velocity  for  a  given  vessel  pressure  and 
projectile  configuration.  The  projectile  velocity  data  from  Table  5-1  is 
indicated  by  the  small  squares  in  the  figure.  The  vertical  bars  about 
each  velocity  data  point  denote  the  accuracy  of  the  measurements  (±4.6 
percent  -  discussed  in  Section  5.2).  The  use  of  carefully  placed 
piezoelectric  pins  proved  to  be  a  reliable  method  of  measuring  the 
projec  .ile  velocity,  with  the  exception  of  Tests  9,  16,  and  17.  The 
TOADS  box  was  faulty  on  Test  9  and  16,  preventing  the  collection  of  data. 
On  Test  17,  the  clock  speed  of  the  TOADS  box  was  inadvertently  set 
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Table  5-1.  Projectile  velocity,  planarity,  and  pressurization  time  from  WES  4-Ft 
Diameter  Vertical  Gas  Gun  tests. 


(Meas.  Vel. )/(Pred.  Vel.)  *  100 
Velocity  based  on  one  set  of  TOA  pins 


PROJECTILE  VELOCITY  (M/SEC) 


to  10  kHz,  rendering  the  data  useless.  Test  20  also  used  piezoelectric 
pins  for  measuring  projectile  velocity.  However,  larger  diameter 
(17.5  mm)  o-rings  were  inadvertently  used  on  the  projectile  for  this 
test.  The  resulting  increased  friction  between  the  projectile  and  the 
barrel  caused  the  projectile  to  hesitate  when  released  by  the  trigger 
mechanism.  This  momentary  hesitation  (approximately  one  second)  allowed 
the  TOADS  clock  to  expire  prior  to  arrival  of  the  projectile  at  the  pin 
locations  at  the  end  of  the  barrel,  thus  prohibiting  a  measurement. 

Tests  3  and  4  used  a  side  way  self- shorting  pin.  The  failure  of  three 
sets  of  pins  out  of  four  triggered  the  change  in  pin  types  to  the 
piezoelectric  pin  now  used. 

The  measured  projectile  velocity  was,  on  the  average,  about  90 
percent  of  the  predicted  value  (see  Table  5-1).  The  discrepancy  between 
the  predicted  and  measured  values  may  perhaps  be  attributed  to  the 
friction  between  the  projectile  o-rings  and  the  barrel  (not  accounted  for 
in  the  mathematical  model) ,  and  the  weight  of  the  water  reaction  mass  at 
the  top  of  the  barrel.  When  the  projectile  was  redesigned,  its  weight 
was  increased  from  761  kg  to  1,450  kg  (1,675  lb  to  3,200  lb).  To  obtain 
the  design  maximum  projectile  velocity,  a  further  increase  of  818  kg 
(1,800  lb)  in  the  reaction  mass  is  required.  The  mass  of  the  water  has 
not  bee\  increased  since  the  redesign  of  the  projectile,  since  that 
would  require  increasing  the  height  of  the  upper  barrel  section  that 
contains  the  water  by  approximately  71  cm  (28  in.). 

The  planarity  of  the  projectile,  as  it  exited  the  barrel,  was 
measured  on  the  gun  performance  tests.  Typically,  the  rotation  was 
between  2  and  5  milliradians .  These  results  are  presented  in  Table  5-1. 

Also  listed  in  Table  5-1  (and  presented  in  Figure  5-5)  is  the  tioa 
required  for  the  air  compressor  to  pressurize  the  gas  gun  to  a  given 
level.  The  air  compressor  currently  in  use  can  only  be  run 
intermittently  at  pressures  greater  than  1,379  kPa  (200  psi),  because  of 
increased  temperatures  experienced  by  the  compressor.  This  causes  a 
slight  increase  in  the  amount  of  time  required  to  pressurize  the  4-ft  gun 
reservoir  to  the  higher  operating  pressures. 
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SECTION  6 


CONCLUSIONS  AND  RECOMMENDATIONS 

6.1  CONCLUSIONS. 

Twenty- two  tests  have  been  conducted  with  the  WES  4-ft  Disaster 
Vertical  Gas  Gun  to  date.  These  tests  were  conducted  over  the  entire 
operating  range  of  the  gun  (0  to  2,070  kPa  (300  psi)).  The  projectile 
velocity  on  these  tests  varied  from  5.4  to  71  a/sec.  The  projectile 
velocities  were,  on  the  average,  about  90  percent  of  the  values  predicted 
by  the  mathematical  model  of  the  gun. 

The  long-range  effects  data  collected  on  these  tests  indicates 
there  is  no  potential  for  damage  to  buildings  in  the  area  from  either 
ground  shock  or  airblast.  The  noise  generated  from  firing  the  gun,  while 
relatively  loud,  has  caused  little  disturbance  in  the  surrounding  area. 

The  redesigned  projectile  survived  all  tests  and  performed  as 
expected.  A  carriage/ impact  plate  configuration,  of  the  type  employed  in 
the  current  design,  should  aid  in  delivering  impact  shock  loadings  with 
varying  pulse  shapes  against  target  specimens  in  future  tests. 

The  consistency  and  predictability  of  the  data  collected  on  these 
tests  will  be  of  significant  benefit  for  predicting  and  evaluating  the 
results  of  future  tests  with  the  gun. 

6.2  RECOMMENDATIONS . 

There  should  be  no  further  requirement  for  the  measurement  of  far- 
field  ground  motion  and  airblast  on  future  tests.  Sound  pressure  level 
measurements  should  be  continued,  as  they  are  simple  and  inexpensive 
measurements  that  will  document  the  actual  long-range  "nuisance"  effect 
levels  of  future  tests,  particularly  since  such  low-level  effects  are 
subject  to  variation  with  daily  weather  conditions. 

Future  tests  will  focus  on  the  stresses  and  particle  velocities 
generated  In  target  samples  of  various  geologic  materials,  which  is  the 
main  purpose  of  the  gun.  Finite  element/finite  difference  calculations, 
to  supplement  the  1-D  elastic  calculations  (not  presented  in  this 
report),  should  be  made  to  aid  in  understanding  the  environments  created 
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in  target  testbeds  beneath  the  gun.  The  calculations  should  characterize 
the  stress  and  velocity  fields  in  a  test  sample.  These  calculations  are 
necessary  to  utilize  the  gun  as  a  dynamic  validation  facility. 
Calculations  to  vary  the  pulse  shape  input  to  a  test  article  are  also 
required  to  fully  utilize  the  facility.  These  calculations  should 
address  axisymetrically  diverging  flows  as  well  as  varying  pulse 
durations . 
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APPENDIX  A 


LIST  OF  ABBREVIATIONS 

A  area  of  barrel 

E  yield  energy  of  TNT  explosive 

KE  kinetic  energy 

F  force  on  projectile  and  reaction  mass 

Lj,  distance  traveled  by  projectile 

L*  distance  traveled  by  reaction  mass 

Mp  projectile  mass 

reaction  mass 

P  peak  pressure  amplitude 

total  pressure  in  system 

P.  effective  pressure 

pinit  initial  pressure  acting  on  projectile  and  reaction  mass 
VImt  reference  pressure 

F.  range  from  source 

V  total  system  volume 

VohM  volume  of  chamber  between  projectile  and  reaction  mass 
VA  peak  vertical  particle  velocity  for  impact 

Vp  projectile  velocity 

V  volume  of  pressure  vessel 

Velp  projectile  velocity 

x  displacement  of  projectile 

y  displacement  of  reaction  mass 

p  density  of  reaction  mass 

7  ratio  of  specific  heats  for  a  gas 
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APPENDIX  B 


STANDARD  OPERATING  PROCEDURE  FOR  THE 
WES  4-FT  DIAMETER  VERTICAL  GAS  GUN 


STANDARD  OPERATING  PROCEDURE  FOR  THE 
WES  4 -FT  DIAMETER  VERTICAL  GAS  GUN 

The  following  standard  operating  procedure  (SOP)  will  be  adhered 
to  when  testing  with  the  WES  4-ft  Diameter  Vertical  Gas  Gun. 

I.  LOADING  THE  PROJECTILE: 

CONTROL  PANEL  STATUS: 

Master  control  switch  on  side  of  control  panel  "ON" 

Poway  direct  switch  on  side  of  control  panel  set  to  "GUN* 

AT  THE  GUN: 

Set  regulator  on  nitrogen  bottle  to  150  psi 
Hook-up  battery  power  to  fire/cock  valve 

Open  manual  valve  Inside  barrel  above  trigger  mechanism  (always 

open) 

Lower  latches  on  trigger  mechanism  by  placing  firo/cock  switch 
in  "FIRE"  position 

Move  projectile  into  place 

Hove  fire/cock  switch  to  the  "COCK"  position 

n. 

Position  fiberglass  sheet  to  support  water  reaction  mass,  add 
foam  board  and  plastic  liner 

Fill  top  of  barrel  with  uater 

Set  regulator  on  air  tank  to  50  psi 

Turn  on  water  to  aftarcoolor 

*****  CLEAR  ALL  PERSONNEL  FROM  SITE  ***** 
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Ill .  REMOTE  CONTROL  FOR  PRESSURIZING  AND  BLEEDING  VESSEL: 
CONTROL  PANEL  STATUS : 


Master  control  switch  on  side  of  control  panel  "ON" 

Power  direct  switch  on  side  of  control  panel  set  to  "CONTROL 

PANEL" 

Master  control  switch  on  front  of  control  panel  set  for  power  to 
left  side  of  control  panel,  left  light  on 

All  val^res  initially  closed 

PRESSURIZE  SEQUENCE: 

•OPEN"  pressure  isolation  valve  in  pressure  line 
•START"  air  compressor 

Monitor  vessel  pressure  with  the  direct-reading  pressure  meter 
(units  of  psi)  mounted  in  control  panel 

"STOP"  compressor  when  desired  pressure  level  is  obtained 

"CLOSE"  pressure  isolation  valve  in  the  pressure  line 

&6BILa^Bgg: 

"OPEN"  bleed  valve  in  pressure  line 

When  vessel  bleeding  is  complete,  "OPEN"  pressure  isolation 
valve  to  bleed  pressure  line  back  to  air  compressor 

"CLOSE"  both  pressure  line  valves 

IV.  REMOTE  CONTROL  FOR  EVACUATING  AND  BLEEDING  BARREL: 

CONTROL  PANEL  STATUS: 

Master  control  switch  on  side  of  control  panel  *QN* 

Power  direct  switch  on  side  of  control  panel  set  to  "CONTROL 

PANEL" 

Master  control  switch  on  front  of  control  panel  set  for  power  to 
left  side  of  control  panel,  loft  light  on 

All  valves  initially  closed 


B-3 


"OPEN"  vacuum  isolation  valve  in  vacuum  line 
“OPEN"  pressure  isolation  valve  in  vacuum  line 
"START"  vacuum  pump 

Monitor  vacuum  in  barrel  with  the  direct -reading  pressure  meter 
(units  of  in.  Hg)  mounted  in  control  panel 

"STOP"  vacuum  pump  when  desired  vacuum  level  is  obtained 

"CLOSE"  vacuum  isolation  valve  in  vacuum  line 

"CLOSE"  pressure  isolation  valve  in  vacuum  line 

BLEED  SEQUENCE: 

"OPEN"  pressure  isolation  valve  in  vacuum  line 
"OPEN"  bleed  valve  in  vacuum  line 

^ii.;n  barrel  bleeding  is  complete,  "OPEN"  vacuum  isolation  valve 
to  bleed  vacuum  line  back  to  vacuum  pump 

"CLOSE"  all  three  vacuum  line  valves 

V.  FIRING  THE  GUN: 

Master  control  switch  on  front  of  control  panel  set  for  power  to 
right  side  of  control  panel,  right  light  on 

Sound  warning  horn  'me  minute  prior  to  firing,  thirty  seconds 
prior  to  firing,  and  three  short  blasts  to  count  down  filial  seconds  prior 
to  firing 


Fire  gun  by  moving  fire/cock  switch  to  "OPEN"  position 

Lower  retaining  collar  of  trigger  mechanism  by  moving  fire/cock 
switch  to  "CLOSED"  position 


VI.  TO  SHUT  DOWN  CONTROL  PANEL: 

Master  control  switch  on  front  of  control  panel  set  for  power  to 
left  side  of  control  panel,  left  light  on 

Bleed  the  pressure  and  vacuum  lines 

"CLOSE"  all  valves 


Move  master  control  switch  on  side  of  control  panel  "OFF1 
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LISTING  OF  FORTRAN  PROGRAM  PARTVEL.FOR 
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C  PROGRAM  NAME:  PARTVEL.FOR 
C  DATE:  APRIL  27,  1991 

C  WRITTEN  BY:  HOWARD  G.  WHITE 
C  USAE  WATERWAYS  EXPERIMENT  STATION 

C  3909  HALLS  FERRY  ROAD 

C  VICKSBURG,  MS  39180-6199 

C  PH:  (601)  634-3391 

C  THIS  PROGRAM  CALCULATES  THE  PEAK  VERTICAL  PARTICLE  VELOCITY 
C  FOR  IMPACTS  AS  A  FUNCTION  OF  RANGE  AND  PROJECTILE  SEISMIC  ENERGY 
C  LEVEL  ASSUMING  A  TNT  EQUIVALENCE  OF  1.41E6  FT- LB/LB- TNT.  THE 
C  PROJECTILE'S  KINETIC  ENERGY  IS  FOUND  FROM  ITS  MASS  AND  VELOCITY. 

C  REFERENCE:  FUNDAMENTAL  EXPERIMENTS  IN  GROUND  SHOCK  PHENOMENOLOGY 
C  WES  MISCELLANEOUS  PAPER  N-73-2 

C  J.G.  WALLACE  AND  J.  FOWLER 

C  MARCH  1973 

C  PARAGRAPH  28,  P36. 

C  VARIABLE  DEFINITION 

C  GC  -  GRAVITATIONAL  CONSTANT  (LBm-FI/LBf-SECA2) 

C  E  -  ENERGY  OF  IMPACT  (LB -TNT) 

C  MP  -  MASS  OF  PROJECTILE  (LSn) 

C  PARTVEL  -  PARTICLE  VELOCITY  AT  RANGE  R  (IN/SEC) 

C  PROJVEL  -  PROJECTILE  VELOCITY  (FT/SEC) 

c  r  -  range  (ft) 

C  TNT  -  TNT  EQUIVALENCE 

C  VHIN  -  MINIMUM  VELOCITY  USED  IN  CALCULATION  (FT/SEC) 

C  VHAX  -  MAXIMUM  VELOCITY  USED  IN  CALCULATION  (FT/SEC) 

INTEGER  VMIN.VMAX 
REAL  HP 

CHARACTER  QUTFILE*15 
10  CONTINUE 

C  ****  DEFINE  CONSTANTS  **** 

GC-32.174 

TNT-1410000.0 

HP-3260. 

VMIN-1 

VHAX-300 

WRITE (* , 50)GC , TNT , HP , VHIN , VHAX 
50  FORMAT (//,5X,' CURRENT  PARAMETER  SETTINGS  ARE  :  ' , 

1//.5X, 'GRAVITATIONAL  CONSTANT  -  ',F6.3,'  (LBta-FT)/(LBf-SEC"2) ' , 
2// ,5X, 'TNT  EQUIVALENCE  -  \F8.0,'  <FT*LB£')/(Lfif-TNT)  ’ , 

3//, 5X, 'PROJECTILE  HASS  -  \F5.0,'  LBo\ 
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4//.5X, 'RANGE  OF  PROJECTILE  VELOCITY  CALCULATED  :  ',12,'  TO  ', 

5  13,'  FT/SEC',//) 

WRITE(*,60) 

60  FORMAT (5X, 'UNITS  OF  THE  RESULTING  PARTICLE  VELOCITY  ARE  IN/SEC') 
WRITE(*,70) 

70  FORMAT ( // , 5X , '  INP-iJT  THE  RANGE  OF  THE  OBSERVATION  POINT  (FT)  :  ') 

READ(*,*)R 

WRITE(*,80) 

80  FORMAT(//,5X, 'OUTPUT  FILE  NAME  :  ') 

READ(*, 90)0UTFILE 
90  F0RMAT(A15) 

OPEN (UNIT-17 , FILE-OUTFILE) 

WRITE(17,*)VMAX 

C  ****  CALCULATION  SECTION  **** 

DO  100  I-VMIN.VMAX 
PROJVEb-I*1.0 

C  ****  EQN  (4)  IN  NOTES  **** 

E-HP*PROJVEL*PROJVEL/(2 . 0*GC*14 10000. ) 

C  ***•>  EQN  (17)  IN  REFERENCE  WITH  UNITS  OF  IN/SEC  **** 

PARTVEL-600 . 0*(E/R**3 )**0 . 5 

WRXTE( 17 , *) PROJVEL , PARTVEL 

100  CONTINUE 

CLOSE  (UNIT*!?) 

350  WRITE(*,*) '  DO  YOU  WISH  TO  RUN  THIS  PROGRAM  AGAIN?* 
URITE(*,*)'  YES  -  1  NO  -  2' 

WRITE (*,*) *  PREFERENCE:  ' 

RCAD(*,*)It>REF 
IF(IPREF.EQ.l)  GOTO  10 

1000  END 
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LISTING  OF  FORTRAN  PROGRAM  GG4PV . FOR 


C  PROGRAM  NAME:  GG4PV . FOR 
C  UPDATED:  OCTOBER  30,  1992 
C  DATE:  JULY  10,  1991 

C  UPDATE  NOTES:  NOW  ALLOW  OPTIONS  OF  USING  A  PRESSURE  PLATE  ON 
C  THE  PROJECTILE  AND  TO  REDUCE  XL  BY  THE  ADDITION  OF  FOAM  AND  AN  IMPACT 
C  PLATE  TO  THE  FRONT  OF  THE  PROJECTILE, 

C  WRITTEN  BY:  HOWARD  G.  WHITE 
C  USAE  WATERWAYS  EXPERIMENT  STATION 

C  3909  HALLS  FERRY  ROAD 

C  VICKSBURG,  MS  39180-6199 

C  PH:  (601)  634-3391 

C  THIS  PROGRAM  CALCULATES  THE  PROJECTILE  VELOCITY  FOR  THE  4- FT 

C  GAS  GUN,  FOR  A  GIVEN  VESSEL  PRESSURE,  OR  RANGE  OF  PRESSURES, 

C  AND  A  GIVEN  PROJECTILE  WEIGHT. 

C  REFERENCE:  AL  OHRT'S  AND  MY  NOTES  ON  “THEORETICAL  ANALYSIS  OF 
C  PRESSURE-VELOCITY  CURVES"  IN  THE  B4-FT  GAS  GUN 

C  CALCULATIONS"  NOTEBOOK 

C  "DESIGN  AND  TESTING  EXPERIENCES  WITH  THE  WES  12 -IN. 

C  AND  4- FT  DIAMETER  VERTICAL  GAS  GUN".  WHITE,  OHRT,  WELCH 

C  AND  JOACHIM.  1991  DNA  INWET  CONFERENCE 

C  VARIABLE  DEFINITION 
C  C  -  CONVENIENT  CONSTANT 
C  COEFF1  -  CONVENIENT  CONSTANT 
C  COEFF2  ~  CONVENIENT  CONSTANT 
C  COEFF3  -  CONVENIENT  CONSTANT 
C  DENRM  -  DENSITY  OF  REACTION  HASS  (SLUGS/FTA3) 

Cl-  DUMMY  VARIABLE 

C  ICALC  -  1  IF  CALC  VEL  AS  FRONT  EXITS  BARREL,  2  IF  REAR 
C  IPPLATE  -  2  IF  PRESSURE  PLATE  IS  USED  ON  PROJECTILE,  1  IF  NOT 
C  OUTFILE  -  OUTPUT  FILENAME  FOR  CALCULATED  DATA 

C  PINC  -  PRESSURE  INCREMENT  FOR  CALCULATIONS  (PSI,  CHANGED  TO  LB/FTA2) 

C  PINIT  -  INITIAL  VESSEL  PRESSURE  (PSI) 

C  PMASSLB  -  PROJECTILE  MASS  (LBn,  CHANGED  TO  SLUGS) 

C  PMASS  -  PROJECTILE  MASS  (SLUGS) 

C  PKAX  -  HAX  PRESSURE  FOR  CALCULATIONS  (PSI) 

C  RATIO  -  RATIO  OF  PROJECTILE  MASS  TO  REACTION  HASS 
C  RBAR  -  RADIUS  OF  BARREL  (FT) 

C  RHASS  -  REACTION  HASS  (SLUGS) 

C  VCHAM  -  VOL  OF  CHAMBER  BETWEEN  PROJECTILE  AND  REACTION  HASS  (FTA3) 

C  VEL  -  VELOCITY  OF  PROJECTILE  (FT/SEC) 

C  VINIT  -  INITIAL  PRESSURIZED  VOL,  BEFORE  HASSES  BEGIN  HOVING  (FTA3) 

C  VTANK  -  VOLUME  OF  TANK  (FTA3) 

C  XFOAM  -  LENGTH  OF  FOAM  ADDED  TO  FRONT  OF  PROJECTILE  (IN  CHANGED  TO  FT) 
C  XL  -  LENGTH  OF  PROJECTILE  TRAVEL  PRIOR  TO  IMPACTING  RATE  PINS  (FT) 

C  XPLATE  -  LENGTH  OF  IMPACT  PLATE  (IN  CHANGED  TO  FT) 

C  YL  -  LENTH  OF  REACTION  HASS  HEIGHT  (AND  TRAVEL)  (FT) 
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CHARACTER  0UTFILE*15 


C  ****  INPUT  PROJECTILE  MASS  **** 

10  WRITE (*,90) 

90  FORMAT <///,5X,  'PROGRAM  GG4PV  (LAST  UPDATED:  30  OCTOBER  1992)') 
WRITE (*,100) 

100  FORMAT (///,5X, 'INPUT  PROJECTILE  MASS  (LBa)  :  ') 

READ(*,*)PMASSLB 

WRITE (*,110) 

110  FORIIAT (/// , 5X , 'INPUT  LENGTH  OF  FOAM  ON  PROJECTILE  (IN)  :  ') 
READ(*,*)XFOAM 

WRITE(*,120) 

120  F0RMAT(///,5X, 'INPUT  LENGTH  OF  IMPACT  PLATE  (IN)  :  ') 

READ(*,*)XPLATE 

125  WRITE(*,130) 

130  FORMAT (///, 5X, 'WAS  A  PRESSURE  PLATE  USED  ON  THE  PROJECTILE?', 

1  /,5X, '1  -  NO  2  -  YES', 

2  // , 5X , 'ENTER  OPTION  :  ') 

READ(*,*)XPPLATE 

IF ( IPPLATE .  LT .  1  .OR.  IPPLATE.GT.2)  THEN 
WRITE(6 ,*) 'YOU  ENTERED  :  '.IPPLATE 
WRITE(6,*) ‘PLEASE  ENTER  A  1  OR  2* 

GOTO  125 
END  IF 

135  WRITE(*,140) 

140  F0RHAT(/// .  5X,  'CALCULATE  THE  VELOCITY  AS : ' , 

1  /,5X,'l  -  THE  FRONT  OF  THE  PROJECTILE  EXITS  THE  BARREL', 

2  /, 5X, '2  -  THE  REAR  OF  THE  PROJECTILE  EXITS  THE  BARREL' , 

3  //,5X, 'ENTER  OPTION  :  ') 

KEAD(*,*) ICALC 

IF( ICALC.LT. I  .OR.  ICALC.GT.2)  THEN 
WRITE (6,*) 'YOU  ENTERED  :  '.ICALC 
WRITE(6,*) 'PLEASE  EN1;ER  A  1  OR  2' 

GOTO  135 
END  IF 

C  ****  INITIAL  CONSTANTS  AND  CALCULATIONS  **** 

PI  -  3.141592654 
RBAR  -  2.0 
DENFLH  -  1.939 
XPOAK-XFOAM/12. 

XPlATE-XPLATE/12 . 

IF  (ICALC  .EQ.  1)  THEN 
XL  -  (128 , 375-XPOAM-XPLATE)/12 . 
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ELSE 

XL  -  166. /12. 

END  IF 

PMASS  -  PMASSLB/32 . 174 
VTANK  -  355.0 
IF  (IPPLATE  .EQ.  1)  THEN 
VCHAM  -  63.7 
ELSE 

VCHAM  -  33.0 
END  IF 

VINIT  -  VTANK+VCHAM 

C  ****  CALCULATE  REACTION  MASS  TRAVEL  DISTANCE  AND  MASS  **** 

ZZ  -  XL*PMASS 
Z  -  P I *RBAR*RBAR*D ENRH 

C  ****  SEE  OHRT’S  NOTES,  P9,  EQN  21,  OR  MY  NOTES,  P17,  EQN  21  **** 
YL  -  SQRT(ZZ/Z) 

RMASS  -  YL*Z 

C  ****  CALCULATE  CONVENIENT  CONSTANTS  **** 

C  ****  SEE  MY  NOTES,  P16,  EQN  12  **** 

RATIO  -  PHASS/RMASS 
C  -  PI*RBAR*RBAR*  ( 1 .  -{-RATIO) 

COEFF1  -  (1 . - (VINIT/(VINIT+XL*C) )**0.4) 

COEFF3  -  ((VTANK/VINIT)**1.4)/O.A*2. 

C  ****  SEE  INUET  PAPER,  P8.  EQN  5  **** 

COEFF2  -  ( COEFF3*VINIT ) / ( PMASS* ( l . + RATIO) ) 

C  ****  CALCULATE  ONE  OR  MANY  POINTS?  **** 

150  URITE(*,160) 

160  FORMAT (///, 5X, •  OPTIONS  * , 

1  //,5X, '1  -  GENERATE  DATA  FILE  FOR  RANGE  OF  PRESSURES*, 

2  /,5X,*2  -  INVESTIGATE  SINGLE  PRESSURE*, 

3  //,5X, ’ENTER  OPTION  :  *) 

READ(*,*)IOFT 

IF(I0PT.LT.1  .OR.  IOPT.GT. 2)  THEN 
URITE(6,*)jYOU  ENTERED  :  \I0PT 
URITE(6,*) 'PLEASE  ENTER  A  1  OR  2* 

GOTO  150 
END  IF 

IF(IOPT.EQ.2)  GOTO  300 

C  ****  SECTION  TO  CALCULATE  FOR  A  RANGE  OF  PRESSURES  **** 

200  WRITER,  210) 

210  FORMAT(///,5X,* INPUT  MAXIMUM  PRESSURE  OF  INTEREST  (PSI)  :  *) 

R£AD(*,*)PHAX 


D-4 


URITE(*,220) 


220  FQRMAT(/,5X, 'INPUT  INCREMENT  FOR  PRESSURE  STEPS  (PSI)  :  ') 
READ(*,*)PINC 

WRITE(*,230) 

230  FORMAT(/,5X, 'ENTER  THE  NAME  OF  THE  OUTPUT  FILE  :  ') 

READ (*, 240) OUTFILE 
240  FORMAT  (A15) 

NPTS  -  INT ( PMAX/PINC ) +1 
PINC  -  PINC*144. 

OPEN (UNIT-10 , FILE-OUTFILE) 

UR1TE(10 , 250)NPTS 
250  FORMAT (15) 

C  ****  BEGIN  LOOP  CALCULATIONS  **** 

DO  260  I-l.NPTS 
PRES  -  (I-1)*PINC 

C  SEE  MY  NOTES,  P16,  EQN  12  **** 

C  ****  SEE  INWET  PAPER,  P8,  EQN  9  **** 

VEL  -  (PRES*C0EFF1*CQEFF2)**9 . 5 
PRES  -  PRES/144. 

WRITE ( 10 , *) PRES , VEL 
260  CONTINUE 

CLOSE(UNIT-IO) 

GOTO  400 

C  **+*  SECTION  TO  CALCULATE  VELOCITY  FOR  A  SINGLE  POINT  **** 

300  WRITE(*,310) 

310  FORMAT(///.5X, 'INPUT  INITIAL  VESSEL  PRESSURE  (PSI)  :  ') 
READ(*,*)PINIT 

C  *+**  SEE  H¥  NOTES,  P16,  EQN  12  **** 

VEL  -  ( PINITH44 .  *COEFFi*COEFF2  )  **0 . 5 
YL  -  YL*12 . 

RMASS  »  RMASS*32.174 
WRITS(* , 320)?INIT, VEL, YL, RHASS , PHASSLB 
320  FORKAT(/,5X, 'INITIAL  PRESSURE  -  ',F4.0,f  PSI', 

1  /, 5X, 'PROJECTILE  VELOCITY  -  ',F4.0,'  FT/SEC', 

2  /,5X, 'REACTION  HASS  HEIGHT  -  \FS.I,'  INCHES', 

3  /, 5X, 'REACTION -MASS  HASS  -  \?S.O,'  LBa\ 

4  /,5X. 'PROJECTILE  MASS  -  \FS.Q,’  LBo') 

C  ****  END  CALCULATIONS  **** 

C  ****  RUN  THE  PROGRAM  AGAIN?  **** 

400  WRITE(*,410) 

410  F08HAT(// , 5X , '  OPTIONS ' , 

1  //,5X,'l  -  RUN  PROGRAM  AGAIN', 

2  /,5X,'2  -  END  PROGRAM', 

3  //,SX, 'ENTER  OPTION  :  ') 
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READ(*,*)I0PT 

IF(IOPT.LT. 1  .OR.  IOPT.GT. 2)  THEN 
WRITE ( 6 , * ) ' YOU  ENTERED  :  '.IOPT 
WRITE(6,*) 'PLEASE  ENTER  A  1  OR  2' 
GOTO  400 
END  IF 

IF  (IOPT. EQ. 1)  GOTO  10 
1000  END 
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APPENDIX  E 


LISTING  OF  FORTRAN  PROGRAM  TOACHECK.FOR 


C  PROGRAM  NAME:  TOACHECK. FOR 
C  DATE:  NOV  25,  1991 


C  WRITTEN  BY: 

C 

C 

C 

C 

C 

C 


HOWARD  G.  WHITE 

USAE  WATERWAYS  EXPERIMENT  STATION 

STRUCTURES  LABORATORY/EXPLOSION  EFFECTS  DIVISION 

ATTN:  CEWES-SE-R 

3909  HALLS  FERRY  ROAD 

VICKSBURG,  MS  391«0-6199 

(601)  634-3391 


C  REFERENCE: 

C 

C 


GAS  GUN  TESTING  NOTEBOOK 
PERFORMANCE  TEST  NO.  5 

NOTES  ON  ANALYSIS  OF  VELOCITY  AND  PLANARITY  MEASUREMENT 


C  GIVEN  THE  SIX  TIME  OF  ARRIVAL  (TOA)  MEASUREMENTS  RECORDED  AT  THE  BASE 
C  OF  THE  BARREL  OF  THE  WES  FOUR  FT  DIAMETER  VERTICAL  GAS  GUN,  THIS 
C  PROGRAM  CALCULATES  THE  VELOCITY  AT  TWO  LOCATIONS  AND  THE  AVERAGE 
C  VELOCITY  OF  THE  TWO.  THE  RELATIVE  DISTANCE  BETWEEN  THE  FOUR  LOCATIONS 
C  ON  THE  BOTTOM  PLATE  OF  THE  PROJECTILE  (USED  AS  A  PLANARITY  CHECK)  ARE 
C  THEN  CALCULATED  BY  FINDING  THE  FIRST  TOA  AND  SUBTRACTING  IT  FROM  THE 
C  THREE  OTHER  TOA  MEASUREMENTS.  THIS  TIME  INCREMENT  IS  THEN  MULTIPLIED 
C  BY  THE  AVERAGE  VELOCITY  TO  DETERMINE  THE  RELATIVE  DISTANCE  BETWEEN  THE 
C  FOUR  POINTS.  THE  XYZ  COORDINATES  ARE  THEN  KNOWN  FOR  EACH  LOCATION. 

C  THREE  OF  THESE  POINTS  ARE  THEN  USED  TO  DETERMINE  THE  EQUATION  OF  THE 
C  PLANE  THAT  PASSES  THRU  THEM.  THIS  EQUATION  IS  THEN  USED  TO  CALCULATE 
C  THE  Z- COORD I NATE  OF  THE  FOURTH  POINT.  THE  PERCENT  DIFFERENCE  BETWEEN 
C  THE  CALCULATED  VALUE  AND  THE  EXPERIMENTAL  VALUE  (INCLUDING  PLACEMENT 
C  TOLERANCES)  IS  THEN  DETERMINED.  THIS  "CHECK1*  IS  HADE  FOR  ALL  FOUR 
C  POINTS  BY  USING  THE  EQUATION  OF  THE  PLANE  THROUGH  THE  OTHER  THREE  PTS. 


C  VARIABLES  USED  IN  THIS  PROGRAM  INCLUDE 


C  COEFF1  **  COEFFICIENT  OF  X  IN  EQN  OF  A  PLANE 

C  COEFF2  -  COEFFICIENT  OF  Y  IN  EQN  OF  A  PLANE 

C  COEFF3  -  COEFFICIENT  OF  Z  IN  EQN  OF  A  PLANE 

C  CONST  -  CONSTANT  TERM  IN  EQN  OF  PLANE 

CD-  DISTANCE  BETWEEN  PIN  LOCATIONS  FOR  VELOCITY  MEASUREMENT  (METERS) 

C  DIA  -  DIAMETER  BETWEEN  TIPS  OF  TOA  PINS  (METERS) 

C  DTMAX  -  MAX  DIFFERENCE  IN  TOA  MEASUREMENTS  (SEC) 

C  DT(*)  -  DIFF  IN  TOA  BETWEEN  A  PIN  LOCATION  AND  TGAMIN  (SEC) 

C  DT*U  -  DIFF  IN  TOA  BETWEEN  A  PIN  LOCATION  AND  TOAMIN  (MICROSECONDS) 

C  IANS  -  ORIENTATION  OF  PINS,  1-20  DEG.  2-HORIZONTAL 
C  IVEL  -  AVERAGE  VELOCITY  EQUALS  EITHER  V8L1  OR  VEL3 
C  PINMAX  -  PIN  LOCATION  OF  HAX  ANGLE  OF  IMPACT.  RELATIVE  TO  PINHIN 
C  PINHIN  -  PIN  LOCATION  OF  TOAMIN 

G  R*-RADIAL  DISTANCE  BETWEEN  LOCATION  OF  TOAMIN  AND  LOCATION  *  (METERS) 
C  RAD  -  RADIUS  BETWEEN  TIPS  OF  TOA  PINS  (METERS) 

C  SM  -  LOCATION  OF  THE  SCRIBE  MARK  FROM  THE  TIP  OF  THE  TOA  PIN  (IN) 

C  TlA  -  TOA  AT  PIN  T-lA  (SEC) 

C  T3A  -  TOA  AT  PIN  T-3A  (SEC) 

C  THETA*-ANGLE  OF  IMPACT  BETWEEN  LOCATION  OF  TOAMIN  AND  LOCATION  *(RAD) 
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C  THETAMAX  -  MAX  ANGLE  OF  IMPACT  (MILLIRADI ANS ) 

C  THETATOL-TOL  ON  ANG  OF  IMPACT  ASSOC  W/  MAX  ANG  OF  IMPACT (MILLIRADIANS ) 
C  TOA(*)  -  TOAs  AT  THE  FOUR  PLANARITY  CHECKPOINTS  (SEC) 

C  TOAMAX  -  MAXIMUM  TOA  VALUE  (SEC) 

C  TOAMIN  -  MINIMUM  TOA  VALUE  (i’EC) 

C  TOL  -  ABSOLUTE  TOLERANCE  ON  PLACEMENT  0?  PI, IS,  (METERS) 

C  TOLMM  -  TOLERANCE  ON  PLACEMENT  Of  PINS,  BEST  WE  CAN  DO  (MILLIMETERS) 

C  TOLPM-PLUS  OR  MINUS  TOLERANCE  ON  PLACEMENT  OF  PINS,  1/2  OF  TOL (METERS) 
C  TOLPMMM-PLUS  OR  MINUS  TOLERANCE  ON  PLACEMENT  OF  PINS,  1/2  OF  TOL  (MM) 

C  VEL1  -  VELOCITY  AT  PIN  LOCATION  1  (M/SEC) 

C  VEL3  -  VELOCITY  AT  PIN  LOCATION  3  (M/SEC) 

C  VELLB-MININUM  VELOCITY  PLUS  ITS  ASSOCIATED  VELOCITY  TOLERANCE  (M/SEC) 

C  VELTOL*  -  TOLERANCE  ON  VELOCITY  AT  PIN  LOCATION  *  (M/SEC) 

C  VELUB-MAXIMUM  VELOCITY  MINUS  ITS  ASSOCIATED  VELOCITY  TOLERANCE  (M/SEC) 
C  VELAVG  -  AVERAGE  VELOCITY  (M/SEC) 

C  VELAVGFPS  -  AVERAGE  VELOCITY  (F/SEC) 

C  X*  -  X  COORDINATE  0*  PLANARITY  CHECK  PINS  (METERS) 

C  X*MM  -  X  COORDINATE  OF  PLANARITY  CHECK  PINS  (MILLIMETERS) 

C  Y*  -  Y  COORDINATE  OF  PLANARITY  CHECK  PINS  (METERS) 

C  Y*MM  -  Y  COORDINATE  OF  PLANARITY  CHECK  PINS  (MILLIMETERS) 

C  Z*  -  Z  COORDINATE  OF  PLANARITY  CHECK  PINS  (METERS) 

C  Z*C  -  CALC  VALUE  OF  Z*  BASES  ON  EQN  OF  PLANE  THRU  OTHER  PINS  (METERS) 

C  Z*CHM  -  CALC  VALUE  OF  Z*  BASED  ON  EQN  OF  PLANE  THRU  OTHER  PINS  (MM) 

C  Z*HAX  -  MAX  VALUE  BETWEEN  1*  AND  (METERS) 

C  Z*MXN  -  MIN  VALUE  BETWEEN  Z*  AND  Z*C  (METERS) 

C  Z*MM  -  Z  COORDINATE  OF  PLANARITY  CHECK  PINS  (MILLIMETERS) 

C  ZDIF*HAX  -  MAX  %  BXFF  BETWEEN  CALC  AND  ACTUAL  VALUE  OP  FIN  LOCATIONS 

C  ZDIF*MIN  -  KIN  %  DIFP  BETWEEN  CAJLC  AND  ACTUAL  VALUE  OF  PIN  LOCATIONS 

C  ZTT>L*  -  «  DIPF  BETWEEN  TOLERANCE  ON  PLACEMENT  &  ACTUAL  PIN  LOCATION. 


DIMENSION  T0A(4),DT(4) 

DOUBLE  PRECISION  TOA, DT, TOAMIN 
CHARACTER  PIRKIE*4 , PINHAX*4 
IVEL  -  0 


C  CONVERT  20  DEGREE  PIN  PLACEMENT  ANGLE  PROS  B&6&KSS  TO  RADIANS 
PI-3 ,141592654 
ANGLE-20.  /180>PI 

WRITE  (*, 5) 

5  FORMAT { ////// ,  $X . ' PROGRAM :  TQACHSGK * . // . SX, *  TEST  NAME :  ' ) 

READ  (*,*) 

WRITE  (*,10) 

10  FORMAT (//.SX, 'INPUT  THE  TOA  DATA  FOR  THE  SIX  PINS', 

2  ///.isxt'«K  msaisa  nu  (sec):  ') 

SRAD(*,*)T1A 
WRITE  {*,20} 

20  FORMAT (//, 15X, 'PIN  NUMBER  T-IB  (SEC) :  *) 

READ  (*,*)T0A<1) 

WRITE  (*,30) 


30 


40 


50 


60 

64 

65 


66 


67 

68 


C 

C 

C 

C 

C 

C 

C 

G 

e 

c 

c 

c 

c 


FORMAT <//,15X, 'PIS  NUMBER  T* 2  (SEC):  ') 

READ  (*,*)T0A(2) 

WRITE  (*,40) 

FORMAT (//,15X, 'PIN  NUMBER  T-3A  (SEC):  *> 

READ  (*,*)T3A 
WRITE  (*,50) 

FORMAT (//,15X, 'PIN  NUMBER  T-3B  (SEC):  ') 

READ  (*,*)TOA(3) 

WRITE  (*,60) 

FORMAT (//,15X, 'PIN  NUMBER  T-4  (SEC):  ') 

READ  (*,*)T0A(4) 

WRITE(*,65) 

FORMAT (///,5X, 'WERE  THE  TOA  PINS  PLACED', 

1  '  AT  A  20  DEGREE  ANGLE  (1)  OR  HORIZONTALLY  (2)  7' , 

2  /,15X, 'ORIENTATION:  ') 

READ (*,*) IANS 

IF  (IANS.SQ.1  .OR.  IANS.EQ.2)  THEN 
GOTO  67 
ELSE 

WRITE  (*,66) 

H}RMAT(/,15X, 'INPUT  A  1  OR  2  PLEASED 
GOTO  64 
ENDIF 

WRITE  (*,68) 

FORMAT (///,SX, 'INPUT  THE  LOCATION  OF  THE  SCRIBE', 

I  /,5X, 'MARK  FROM  THE  TIP  OF  THE  TOA  PIN  (IN):  ') 

READ(*,*)SH 

AN  ASSUMPTION  HERS  IS  THAT  THE  X  AMO  Y  LOCATIONS  OF  THE  PINS  DOES 
NOT  CHANGE.  IE,  THE  CENTER  OF  THE  PROJECTILE  IS  IN  THE  CENTER  OF  THE 
BARREL  AT  THE  TIME  OF  PIN  IMPACT. 


TOLERANCES  014  THE  VELOCITY'  AND  PLANARITY  MEASUREMENTS 


FACTOR 


TOLERANCE  (+0*1  -) 


20  DECREE  ANGLE 


PIN  PLACEMENT  0.005" 
RECEPTACLE  LOCATION  0.001" 
PIN  HOLES  IN  RECEPTACLE  0.002" 
FLATNESS  OF  IMPACT  PLATE  0.005" 
PROJECTILE  ASSEMBLY  0.010“ 
TOTAL  +  OR  .  0.023° 
ABSOLUTE  TOLERANCE  0.046" 


+  OR 


HORIZONTAL 

0.000" 

0.G01" 

0.000" 

0.005" 

0.010" 

*  0.016* 
0.032* 


C  DIA  IS  DIAMETER  BETWEEN  PINS  RESULTING  FROM  A  PLACEMENT  MARK  AT  SM" . 
C  THE  ANGLE  AT  WHICH  THE  PIN  IS  PLACED  IS  20  DEGREES. 


IF(IANS.EQ.2)  THEN 
TOL-O.Q32 
DIA-48.  -  (2.*SK) 

ELSE 

TOL-O.046 

DIA-48.  -  (2 .*(SH*COS(ANGLE) ) ) 


ENDIF 


C  CONVERT  TOLERANCE  TO  MILLIMETERS  AND  METERS  AND  DIAMETER  TO  METERS 
TOLMM-TOL*25.4 
TQL-TQL*0 . 0254 
TOLPM-TGL/2 . 

TOLPMMM-TOLMM/2 . 

DIA-DIA/12./3.28 
RAD-DIA/2 . 

C  CALC  THE  VELOCITY  AT  LOCATIONS  1  AND  3 
C  THE  DISTANCE  BETWEEN  THE  PINS  IS  1",  CHANGED  TO  METERS 

D  -  1./12./3.28 
VEL1  -  D/(T0A{1) -T1A) 

VELTOL1  -  TOL/ (TOA ( 1 ) - TIA) 

VEL3  -  D/(T0A(3) -T3A) 

VELTOL3  -  TOL/(TOA(3) -T3A) 

WRITE  (*,70)  VEL1 , VELTOLl , VEL3 , VELTOL3 
70  FORMAT (///, 5X, 'THE  VELOCITY  AT  PIN  LOCATION  1  IS 

1  /, 5X, '  PLUS  OR  MINUS 

2  //, 5X, 'THE  VELOCITY  AT  PIN  LOCATION  3  IS 

3  / , 5X , '  PLUS  OR  MINUS 


' , F6 . 2 , '  M/SEC', 
\F6.3,'  M/SEC', 
'  ,F6.2,'  M/SEC', 
' , F6 . 3 , '  M/SEC') 


C  CHECK  TO  SEE  IF  THE  AVG  VELOCITY  WITHIN  THE  BOUNDS  OF  THE  TOLERANCES 


VELAVG  -  (VEL1+VEL3 ) /2 . 
VELMIN-AMIN1 ( VEL1 , VEL3 ) 
VELMAX-AMAX1 ( VELl , VEL3 ) 
IF  ( VELMIN . EQ , VELl )  THEN 
VELLB-VELMIN+VELTOL.1 
VELUB-VELMAX-VELTOL3 
ELSE 

VELLB-VELMIN+VELTOL3 
VELUB-VELKAX  -  VELTOLl 
ENDIF 


VELAVGFPS  -  VELAVG*3 . 28 


72 


IF( VELAVG . LT . VELLB  .OR.  VELAVG. CT.VELUB)  THEN 
WRITE (* , 72 ) VELAVG , VELAVGFPS , VELUB ,VELLB 
FORMAT ( // , 5X , ' THE  AVERAGE  VELOCITY:  ',F6.2,'  M/SEC', 

L'  ',F5.I,'  FPS'  ,/,5X,  'IS  WITHIN  THE  UPPER  BOUND:  '.F6.2, 


2'  M/SEC' ,/,'X, 'AND  THE  LOWER  BOUND: 


' , F6 . 2 , '  M/SEC', 


3  /,5X, 'OF  THE  VELOCITY  MEASUREMENTS. ' 

4  /, 5X, 'NOTE:  THE  UPPER  BOUND  IS  THE  MAXIMUM  VELOCITY  MINUS' 

5  / , 5X , '  ITS  ASSOCIATED  TOLERANCE. ' 

6  / , 5X , ' NOTE :  THE  LOWER  BOUND  IS  THE  MINIMUM  VELOCITY  PLUS', 

7  / , 5X , '  ITS  ASSOCIATED  TOLERANCE. ') 


ELSE 

WRITE (* , 74 ) VELAVG , VELAVGFPS , VELUB , VELLB 
FORMAT ( // , 5X , 'THE  AVERAGE  VELOCITY: 
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' , F6 . 2 , '  M/SEC', 


1' 

2' 

3 

4 

5 

6 
7 


' , F5 . 1 , '  FPS' ,/,5X, 'IS  NOT  WITHIN  THE  UPPER  BOUND:  \F6.2, 
M/SEC ' , / , 5X, ' AND  THE  LOWER  BOUND :  ' , F6 . 2 , '  M/SEC ' , 

/,5X, 'OF  THE  VELOCITY  MEASUREMENTS. ' 

/,5X, 'NOTE:  THE  UPPER  BOUND  IS  THE  MAXIMUM  VELOCITY  MINUS', 
/ , 5X, '  ITS  ASSOCIATED  TOLERANCE. ' 

/,5X, 'NOTE:  THE  LOWER  BOUND  IS  THE  MINIMUM  VELOCITY  PLUS', 

/ , 5X , '  ITS  ASSOCIATED  TOLERANCE. ') 


75  WRITE (*, 76 )VELl,VEL3 

76  FORMAT (.//,5X, 'CONTINUE  THE  PROGRAM  WITH  EITHER: ',//, 15X, 
l'VELl  -  ' ,F6.2,/,8X, 'OR' ,5X, 'VEL3  -  ' ,F6.2,//,5X, 'AS  THE  AVERAGE', 
2'  VELOCITY. ‘ ,//,5X, 'ENTER  A  "1"  OR  A  "3”:  ') 

READ(*,*)IVEL 

IF(IVEL.EQ.l  .OR.  IVEL.EQ.3)  THEN 
GOTO  79 
ELSE 

WRITE(*,78) 

78  FORMAT(//,5X, 'INPUT  A  1  OR  3  PLEASE!') 

GOTO  75 

ENDIF 

79  ENDIF 


IF(IVEL.EQ.l)  THEN 
VELAVG  -  VELl 
ELS  E I F ( I V  EL . EQ . 3 ) THEN 
VELAVG  -  VEL3 
ENDIF 

WRITE(*,80) 

80  FORMAT(//,5X, 'PRESS  ANY  KEY  TO  CONTINUE',//) 
READ(*,*) 


C  CALCULATE  THE  DIFFERENCE  IN  TOA  (DT(*) )  AND  DEFINE  Z*  BASED  ON  VELAVG 
C  ZTOL*  IS  MADE  A  LARGE  NUMBER  WHEN  THE  PERCENT  DIFFERENCE  IS  CHECKED 
C  AGAINST  THE  LOCATION  WHERE  TOAMIN  HAS  BEEN  DETERMINED.  THIS  LOCATION 
C  HAS  A  Z  VALUE  EQUAL  TO  ZERO.  THE  PERCENT  DIFFERENCE  CANNOT  BE 
C  DETERMINED  WHEN  COMPARING  A  CALCULATION  TO  ZERO. 

C  THE  LARGE  NUMBER  WILL  PRINT  TO  THE  SCREEN  AS  ****** 


C  FIND  THE  MINIMUM  AND  MAXIMUM  TOA  VALUE 
TOAMIN  -  TOA(l) 

DO  90  1-2,4 

TOAMIN  -  DMIN1 (TOAMIN , TOA( I ) ) 

90  CONTINUE 

IF  (TOAMIN. EQ.TOA(l))  THEN 
PINMIN-'T-lB' 

DT(I)~Q. 

DT(2)-TOA(2) -TOA(l) 

DT(3)-TOA(3) -T0A(1) 

DT(4)-T0A(4) -TOA(l) 

Xl-O. 
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Yl-O. 

Zl-O. 

X2— RAD 
Y2-RAD 

Z2-VELAVG*DT ( 2 ) 

X3-0. 

Y3-D1A 

Z3-VELAVG*DT ( 3 ) 

X4-RAD 

Y4-RAD 

Z4-VELAVG*DT ( 4 ) 

ELSE  IF  (TOAMIN . EQ . TOA ( 2) )  THEN 
PINMIN-'T-2' 

DT(l)-TOA(l) -TOA(2) 

DT(2)-0. 

DT(3)-TOA(3)-TOA(2) 

DT(4)-T0A(4)-T0A<2) 

Xl— RAD 
Yl-RAD 

Z1-VELAVG*DT(1) 

X2-0. 

Y2-Q. 

Z2-0. 

X3— RAD 
Y3-RAD 

Z3«VELAVG*DT ( 3  ) 

X4-0. 

Y4-DIA 

Z4-VELAVG*DT<4) 

ELSE  IF  (TOAMIN. EQ.TOA( 3))  THEN 
PINMIN-'T-3B' 

DT(l)-TOA(I) -TOA(3) 
DT(2>-TOA(2)-TOA(3) 

DT(3)-0. 

DT(4)-TOA(4)-TOA(3) 

Xl-0. 

Yl-DIA 

Z1-VELAVG*DT ( I ) 

X2-RAD 

Y2-RAD 

Z2-VELAVG*DT ( 2 ) 

X3-0. 

Y3-0. 

Z3-0. 

X4— RAD 
Y4-RAD 

Z4-VELAVG*DT<4) 

ELSE 

PINMIN“'T-4' 

DT(1)-TOA(1)-TOA(4) 

DT(2)-TOA(2)-TOA(4) 

uT(3)-TOA(3)-TOA(4) 
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DT(4)-0. 

XI— RAD 
Yl-RAD 

Z1-VELAVG*DT ( 1 ) 
X2-0. 

Y2-DIA 

22-VELAVG*DT(2) 

X3-RAD 

Y3-RAD 

Z3-VELAVG*DT(3) 

X4-0. 

Y4-0. 

Z4-0. 

ENDIF 


C  CONVERT  DT(*)  TO  MICROSECONDS  TO  PRINT  TO  SCREEN 
DT1U  -  DT(I)  *  1000000. 

DT2U  -  DT(2)  *  1000000. 

DT3U  -  DT(3)  *  1000000. 

DT4U  -  DT(4)  *  1000000. 

C  CONVERT  XYZ  COORDINATES  TO  MILLIMETERS  TO  PRINT  TO  SCREEN 
X1MM  -  Xl*1000. 

X2MM  -  X2*1000 . 

X3MM  -  X3*1000. 

X4MM  -  X4*1000 . 

Y1MM  -  Yl*1000. 

Y2MM  -  Y2*1000 . 

Y3MM  -  Y3*1000. 

Y4MM  -  Y4*1000 , 

Z1MM  -  Z1*1GG0. 

Z2MM  -  Z2*1000 . 

Z3MM  -  Z3*1000 . 

Z4MM  -  Z4*1000. 


WRITE (* , 93 ) DT1U , DT2U , DT3U , DT4U 

93  FORMAT(/,5X, 'THE  DIFFERENCE  IN  TOA  FOR  THE  PLANARITY  CHECK' 

1  '  PINS  IS', 

2  //,15X, 'PIN  T-1B:  ' ,F6.1, '  MICRO -SECONDS' , 

3  /,15X, 'PIN  T-2:  ' ,F6.1, '  MICRO- SECONDS ' , 

4  /.15X, 'PIN  T-3B:  ' ,F6 .1, '  MICRO -SECONDS' , 

5  /,15X, 'PIN  T-4:  ' ,F6.1, '  MICRO- SECONDS') 


WRITE  ( * ,  95 )  X1MM ,  Y1MM ,  Z1MM.X2MM,  Y2MM ,  Z2MM ,  X3MM ,  Y3MM ,  Z3MM , 

1  X4MM ,  Y4MM ,  7.4MM ,  VELAVG ,  TCLPMMM 

95  FORMAT ( // , 5X , ' THE  XYZ  COORDINATES  OF  THE  IMPACT  PLATE  ARE:' 

1  //,21X,'X',10X,'Y',8X,'Z', 

2  /, 5X, 'PIN  T-1B:  ',F8.3,3X,F8.3,1X,F8.3,' 

3  /,5X, 'PIN  T-2:  ' ,F8.3,3X,F8.3,1X,F8.3, ' 

4  /,5X, 'PIN  T-3B:  ' .FS.S.SX.FS.S.IX.FS.S, ' 

5  /,5X, 'PIN  T-4:  '.FS.S.SX.FS.S.lX.FS.S,' 

6  //,5X, 'THE  Z  COORDINATES  OF  THE  IMPACT  PLATE  ARE  BASED  ON' 

7  / , 5X , ' THE  AVERAGE  VELOCITY  ' , F6 . 2 , '  M/SEC ' , 


MILLIMETERS', 

MILLIMETERS', 

MILLIMETERS', 

MILLIMETERS', 
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8  //, 5X, 'THE  TOLERANCE  ON  PLACEMENT  OF  THE  TOA  PINS  IS', 

9  /,5X, 'PLUS  OR  MINUS  '.F5.3,'  MILLIMETERS') 

WRITE (*,80) 

READ(*,*) 

C  PERFORM  CALCS  TO  DETERMINE  MAX  ANGLE  OF  ROTATION  OF  IMPACT  PLATE 
C  ALSO  DETERMINE  TOLERANCE  ON  ANGLE  OF  ROTATION 

IF  ( TOAMIN . EQ . TOA ( 1 ) )  THEN 
R2-SQRT (X2*X2  +  Y2*Y2) 

THETA2-ATAN ( Z2/R2 ) 

R3-SQRT (X3*X3  +  Y3*Y3) 

THETA3— ATAN ( Z3/R3 ) 

R4-SQRT (X4*X4  +  Y4*Y4) 

THETA4-ATAN { Z4/R4 ) 

THETAMAX-AMAX1 (THETA2 , THETA3 , THETA4) 

I F ( THETAMAX . EQ . THETA2 )  THEN 
THETATOL-ATAN ( TOL/R2 ) 

ELSE  IF  (THETAMAX. EQ.THETA3)  THEN 
THETATOL-ATAN (TOL/R3 ) 

ELSE 

THETATOL-ATAN ( T0L/R4 ) 

ENDIF 

ELSE  IF  (TOAMIN. EQ.TOA(2))  THEN 
Rl-SQRT (X1*X1  +  Y1*Y1) 

THETAl-ATAN (Zl/Rl ) 

R3-SQRT (X3*X3  +  Y3*Y3) 

THETA3-ATAN ( Z3/R3 ) 

R4-SQRT (X4*X4  +  Y4*Y4) 

THETA4-ATAN ( Z4/R4 ) 

THETAMAX-AMAX1 (THETAl , THETA3 , THETA4) 

I F ( THETAMAX . EQ . THETA 1 )  THEN 
THETATOL-ATAN ( TOL/R1 ) 

ELSE  IF  (THETAMAX. EQ.THETA3)  THEN 
THETATOL-ATAN ( TOL/R3 ) 

ELSE 

THETATOL-ATAN ( TOL/R4 ) 

ENDIF 

ELSE  IF  (TOAMIN . EQ . TOA ( 3 ) )  THEN 
Rl-SQRT (X1*X1  +  Y1*Y1) 

THETAl-ATAN (Zl/Rl) 

R2-SQRT(X2*X2  +  Y2*Y2) 

THETA2-ATAN ( Z2/R  2 ) 

R4-SQRT (X4*X4  +  Y4*Y4) 

THETA4-ATAN ( Z4/R4 ) 

THETAMAX-AMAX1 (THETAl , THETA2 , THETA4) 

IF(THETAMAX . EQ . THETAl )  THEN 
THETATOL-ATAN ( TOL/Rl ) 

ELSE  IF  (THETAMAX. EQ.THETA2)  THEN 
THETATOL-ATAN ( TOL/R2 ) 

ELSE 
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THETATOL-ATAN ( T0L/R4 ) 

END  IF 
ELSE 

Rl-SQRT (X1*X1  +  Y1*Y1) 

THETA1-ATAN ( Zl/Rl ) 

R2-SQRT (X2*X2  +  Y2*Y2) 

THETA2-ATAN ( Z2/R2 ) 

R3-SQRT (X3*X3  +  Y3*Y3) 

THETA3-ATAN ( Z3/R3 ) 

THETAMAX-AMAX1 (THETAl , THETA2 , THETA3 ) 

IF ( THETAMAX . EQ . THETAl )  THEN 
THETATOL-ATAN (TOL/R1 ) 

ELSE  IF  (THETAMAX. EQ.THETA2)  THEN 
THETATOL-ATAN (TOL/R2 ) 

ELSE 

THETATOL-ATAN ( TOL/R3 ) 

ENDIF 

ENDIF 

IF(THETAMAX.EQ. THETAl)  THEN 
PINMAX-'T-1B' 

ELSE  IF  (THETAMAX. EQ.THETA2)  THEN 
PINMAX-'T-2' 

ELSE  IF  (THETAMAX. EQ.THETA3)  THEN 
PINMAX-'T*3B' 

ELSE 

PINMAX-'T-4* 

ENDIF 

C  CHANGE  THETAMAX  AND  THETATOL  FROM  RADIANS  TO  DEGREES 
THETAMAX  -  THETAMAX*1000 . 

THETATOL  -  THETATOL*10Q0 . 

WRITE(* , 98) THETAMAX , THETATOL, PINMIN , PINMAX 
98  FORMAT (/,2X,’ THE  MAX  ANGLE  OF  ROTATION  OF  THE  IMPACT  PLATE  IS:  ' 

1  F6 . 3 ,  '  MILLIRADIANS' , 

2  /,2X, 'THE  TOLERANCE  ON  THIS  ANGLE  OF  ROTATION  IS  PLUS', 

3  '  OR  MINUS:  ' ,F6.3, '  MILLIRADIANS', 

4  /,2X, 'THIS  ANGLE  OCCURS  BETWEEN  PINS  \A4,'  AND  ' ,A4) 

WRITE(*,80) 

READ(*,*) 

C  PERFORM  CALCS  TO  DETERMINE  CONSISTENCY  OF  TOA  MEASUREMENTS. 

C  THE  Z  COORDINATES  ARE  CALCULATED  BY  ASSUMING  A  PROJECTILE  VELOCITY 
C  EQUAL  TO  VELAVG.  THE  THREE  BASE  POINTS  ARE  ASSUMED  TO  BE  CORRECT. 

C  THE  TOLERANCE  IS  CONSIDERED  ONLY  FOR  THE  POINT  BEING  CHECKED.  THIS 
C  TOLERANCE  IS  NOT  THE  ABSOLUTE  TOLERANCE  (TOL)  BUT  ONE-HALF  THE 
C  ABSOLUTE,  IE,  THE  PLUS  OR  MINUS  TOLERANCE  (TOLPM) . 

C  CHECK  PIN  1  BY  USING  EQN  OF  PLANE  THRU  PINS  2,  3,  &  4 
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100 


C0EFF1  -  (Y2-Y4)*(Z3-Z4)  -  (Y3-Y4)*(Z2-Z4) 

C0EFF2  -  (X2-X4)*(Z3-Z4)  -  (X3-X4)*(Z2-Z4) 

C0EFF3  -  (X2-X4)*(Y3-Y4)  -  (X3-X4)*(Y2-Y4) 

CONST  -  - (X4*C0EFF1)  +  (Y4*C0EFF2)  -  (Z4*COEFF3) 

Z1C  -  ( (Y1*C0EFF2)  -  (Xl*CGEFFi)  -  C0NST)/C0EFF3 
Z1CHM  -  Z1C  *  1000. 

WRITE (* , lOOJCOEFFl , C0EFF2 , C0EFF3 , CONST 

FORMAT (//,18X, '****  CHECK  LOCATION  OF  TOA  PIN  T-1B  ****', 

1  //,5X,'THE  EQUATION  OF  THE  PLANE  THROUGH  TOA  PINS  2,  3, 

2  '  &  4  IS:', 

3  //.10X.E11.5,'  X  +  '.Ell. 5,'  Y  +  \E11.5,'  Z  +  '.EU.5, 


-  O') 


ZIMIN  -  AMINl(Zl.ZlC) 

IF (ZIMIN. EQ. Zl)  THEN 

C  CALCULATED  VALUE  OF  Zl  IS  GREATER  THAN  THE  EXPERIMENTAL  VALUE  OF  Zl 
ZDIF1MIN  -  ( (Z1C/(Z1+T0LPM) )  -  1.)  *  100. 

IF(ZDIF1MIN.LE.0.)  THEN 
WRITE (*,102) Z1CMM , TOLPMMM , Z1MH 

102  F0RMAT(//,5X, 'THE  CALCULATED  VALUE  OF  PIN  T- IB:  '.F6.3, 

1  '  MILLIMETERS', 

2  /,5X, 'IS  WITHIN  THE  PLACEMENT  TOLERANCE:  '.F5.3, 

3  '  HILLIMETERS ' , 

4  /,5X, 'OF  THE  EXPERIMENTAL  VALUE:  ',F6.3, 

5  '  MILLIMETERS ' ) 

ELSE 

ZDIF1MAX  -  ABS((Z1C/(Z1-T0LPM))  -  1.)  *  100. 

WRITE(* , 104) Z1CMM , TOLPMMM , Z1MH , ZDIF1MIN , ZDIFIMAX 
104  FORMAT (//,5X, 'THE  CALCULATED  VALUE  OF  PIN  T-1B:  \F6.3, 

1  '  MILLIMETERS', 

2  /,5X,'IS  NOT  WITHIN  THE  PLACEMENT  TOLERANCE:  \FS.3, 

3  '  MILLIMETERS', 

4  /,5X, 'OF  THE  EXPERIMENTAL  VALUE:  ',F6.3, 

5  '  MILLIMETERS', 

6  // , 5X , 1  THE  MINIMUM  PERCENT  DIFFERENCE  BETWEEN  THE  CALCULATED', 

7  /, 5X, 'AND  EXPERIMENTAL  VALUES  (CONSIDERING  THE  PLACEMENT', 

8  /,5X, 'TOLERANCE)  IS:  ',F7.2,'  PERCENT', 

9  // , 5X , ' THE  MAXIMUM  PERCENT  DIFFERENCE  BETWEEN  THE  CALCULATED' , 

1  /,5X, 'AND  EXPERIMENTAL  VALUES  (CONSIDERING  THE  PLACEMENT', 

2  /,5X, 'TOLERANCE)  IS:  ',F7.2,'  PERCENT') 

ENDIF 

C  CALCULATED  VALUE  OF  Zl  IS  LESS  THAN  THE  EXPERIMENTAL  VALUE  OF  Zl 
C  THE  CALCULATED  VALUE  OF  Zl  IS  GREATER  THAN  ZERO 


106 


ELSE  IF  (Z1C.GT.0)  THEN 
ZDIF1MIN«(((Z1-T0LPM)/Z1C)  -  1.)  *  100. 

IF(ZDIFIMIN.LE.O)  THEN 
WRITE(*,106)Z1CMH, TOLPMMM, Z1MH 

FORMAT (//,5X, 'THE  CALCULATED  VALUE  OF  PIN  T-1B:  ',F6.3, 
1  '  HILLIMETERS * , 


E-U 


2 

3 

4 

5 


/,5X, 'IS  WITHIN  THE  PLACEMENT  TOLERANCE:  ',F5.3, 
'  MILLIMETERS', 

/, 5X, 'OF  THE  EXPERIMENTAL  VALUE:  '.F6.3, 

'  MILLIMETERS') 


ELSE 

ZDIFlMAX-( ( (Zl+TOLPM) /Z1C)  -  1.)  *  100. 

WRITE ( * , 108 ) Z1CMM , TOLPMMM , Z1MM , ZDIF1MIN , ZDIF1MAX 
108  FORMAT(//,5X, 'THE  CALCULATED  VALUE  OF  PIN  T- IB:  '.F6.3, 

1  '  MILLIMETERS', 

2  /,5X, 'IS  NOT  WITHIN  THE  PLACEMENT  TOLERANCE:  \F5.3, 

3  '  MILLIMETERS’, 

4  /,5X, 'OF  THE  EXPERIMENTAL  VALUE:  ',F6.3, 

5  '  MILLIMETERS', 

6  //, 5X, 'THE  MINIMUM  PERCENT  DIFFERENCE  BETWEEN  THE  CALCULATED' 

7  /,5X, 'AND  EXPERIMENTAL  VALUES  (CONSIDERING  THE  PLACEMENT', 

8  /,5X, 'TOLERANCE)  IS:  \F7.2,'  FERCENT', 

9  // , 5X , ' THE  MAXIMUM  PERCENT  DIFFERENCE  BETWEEN  THE  CALCULATED’ 

1  /, 5X, 'AND  EXPERIMENTAL  VALUES  (CONSIDERING  THE  PLACEMENT', 

2  /,5X, 'TOLERANCE)  IS:  ',F7.2,'  PERCENT') 

ENDIF 


ELSE  IF( (Z1C- (Z1-T0LPM) ) ,GT.0)THEN 
WRITE (* , 200) Z1CMM , TOLPMMM , Z1MM 

200  FORMAT (// , 5X , ' THE  CALCULATED  VALUE  OF  PIN  T- IB:  '.F6.3, 

1  '  MILLIMETERS ' , 

2  /, 5X, 'IS  WITHIN  THE  PLACEMENT  TOLERANCE:  \F5.3, 

3  '  MILLIMETERS', 

4  /.5X, 'OF  THE  EXPERIMENTAL  VALUE:  '.F6.3, 

5  '  MILLIMETERS') 

ELSE  IF(ABS(ZlC).LE.ABS(Zi«TOLPM))THEN 
ZDIF1MIN  -AES(1  *  ((Zi-TOLPK)/ZlC))  *  100 
ZDIF1MAX  -ABS(1  -  ( (Z1+TQLPM)/Z1C))  *  100 
WRITE (* , 204) Z1CMM , TOLPMMM, Z1MM , ZDIF1KIN , ZDIF1MAX 
ELSE 

ZDIF1MIN  -ABS ( (Z1C/(Z1 -TOLPM) )  -  1)  *  100 

IF(ABS(Z1C).LE.ABS(21+T0LPM))THEN 
ZDIF1MAX  -ABS(1  -  ((Z1+T0LPM)/Z1C))  *  100 
ELSE 

ZDIF1MAX  -ABS((Z1C/(Z1+T0LPM))  *  1)  *  100 
ENDIF 


WRITE (* , 204) Z1CMM , TOLPMMM , Z1MM , ZDIF1HIN , ZDIFIMAX 
204  FORNAT(//,5X, 'THE  CALCULATED  VALUE  OF  PIN  T-1B:  \F6.3, 

1  '  MILLIMETERS', 

2  /,5X, 'IS  NOT  WITHIN  THE  PLACEMENT  TOLERANCE:  ',F5.3, 

3  ’  MILLIMETERS’, 

4  /,5X,' OF  THE  EXPERIMENTAL  VALUE:  ',F6.3, 

5  '  MILLIMETERS’, 

6  // , 5X, 'THE  MINIMUM  PERCENT  DIFFERENCE  BETWEEN  THE  CALCULATED' 
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7  /,5X, 'AND  EXPERIMENTAL  VALUES  (CONSIDERING  THE  PLACEMENT', 

8  /,5X, 'TOLERANCE)  IS:  \F7.2,'  PERCENT', 

9  //, 5X, 'THE  MAXIMUM  PERCENT  DIFFERENCE  BETWEEN  THE  CALCULATED', 

1  /, 5X, 'AND  EXPERIMENTAL  VALUES  (CONSIDERING  THE  PLACEMENT', 

2  /,bX, 'TOLERANCE)  IS:  '.F7.2,'  PERCENT') 


END  IF 

WRITE (*,80) 

READ(*,*) 

C  CHECK  PIN  2  BY  USING  EQN  OF  PLANE  THRU  PINS  1,  3,  64 

COEFF1  -  (Y4-Y1)*(Z3-Z1)  -  (Y3-Yl)*(Z4-Zl) 

COEFF2  -  (X4-X1)*(Z3-Z1)  -  (X3-X1)*(Z4*Z1) 

COEFF3  -  (X4-X1)*(Y3-Yl)  -  (X3-XI)*(Y4-Y1) 

CONST  -  - (X1*C0EFF1)  +  (Y1*C0EFF2)  -  (Zl*COEFF3) 

Z2C  -  ( (Y2*COEFF2)  -  (X2*C0EFF1)  -  CONST)/COEFF3 
Z2CMM  -  Z2C  *  1000. 

WRITE (* , 110 ) COEFFl , COEFF2 , COEFF3 , CONST 
110  FORMAT ( // , 18X , ' ****  CHECK  LOCATION  OF  TOA  PIN  T-2  ****', 

1  //, 5X, 'THE  EQUATION  OF  THE  PLANE  THROUGH  TOA  PINS  1,  3,', 

2  '  &  4  IS: ' , 

3  //,10X,E11.5, '  X  +  r  ,E11.5, '  Y  +  '.E11.5,'  Z  +  \E11.5,'  «  O') 
Z2MIN  «  AMIN1(Z2,Z2C) 

IF(Z2MIN . EQ.Z2)  THEN 

C  CALCULATED  VALUE  OF  Z2  IS  GREATER  THAN  THE  EXPERIMENTAL  VALUE  OF  Z2 
ZDIF2HIN  -  ( (Z2C/(Z2+TOLPM) )  -  1.)  *  100. 

IF(ZDIF2HIN.LE.O. )  THEN 
WRITE ( * , 1 1 2 ) Z2CMM . TOLPMMM , Z2HH 

112  FORMAT(//,5X, 'THE  CALCULATED  VALUE  OF  PIN  T-2:  ',F6.3, 

1  '  MILLIMETERS', 

2  /, 5X, 'IS  WITHIN  THE  PLACEMENT  TOLERANCE:  *,F5.3, 

3  '  MILLIMETERS', 

4  /,5X, 'OF  THE  EXPERIMENTAL  VALUE:  \F6.3, 

5  '  HILLIMETERS' ) 

ELSE 

ZDIF2HAX  -  ABS ( (Z2C/(Z2 -TOLPH) )  -  1.)  *  100. 

WRITE(* , 114)Z2CMM , TOLPMMM , Z2MM , ZDIF2NIN , ZDIF2MAX 
114  F0RMAT(//,5X, 'THE  CALCULATED  VALUE  OF  PIN  T-2:  \F6.3, 

1  '  MILLIMETERS', 

2  /, 5X, ' IS  NOT  WITHIN  THE  PLACEMENT  TOLERANCE:  \F5.3, 

3  •  HILLIMETERS', 

4  /,5X, 'OF  THE  EXPERIMENTAL  VALUE:  \F6.3, 

5  '  HILLIMETERS', 

6  //,5X, 'THE  MINIMUM  PERCENT  DIFFERENCE  BETWEEN  THE  CALCULATED', 

7  /.SX.'AND  EXPERIMENTAL  VALUES  (CONSIDERING  THE  PLACEMENT', 

8  /,5X, 'TOLERANCE)  IS:  \F7.2,'  PERCENT', 
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9  //, 5X, 'THE  MAXIMUM  PERCENT  DIFFERENCE  BETWEEN  THE  CALCULATED', 

1  / , 5X , ' AND  EXPERIMENTAL  VALUES  (CONSIDERING  THE  PLACEMENT', 

2  /,5X, 'TOLERANCE)  IS:  ',F7.2,'  PERCENT') 

ENDIF 

C  CALCULATED  VALUE  OF  Z2  IS  LESS  THAN  THE  EXPERIMENTAL  VALUE  OF  Z2 

ELSE  IF  (Z2C.GT.0)  THEN 

ZDIF2MIN- ( ( ( Z2 - TOLPM) /Z2C )  -  1.)  *  100. 

IF(ZDIF2MIN.LE.O)  THEN 
WRITE (*, 210) Z2CMM , TOLPMMM , Z2MM 

F0RMAT(//,5X, 'THE  CALCULATED  VALUE  OF  PIN  T- IB:  \F6.3, 

'  MILLIMETERS', 

/,5X, 'IS  WITHIN  THF.  PLACEMENT  TOLERANCE:  ',F5.3, 

'  MILLIMETERS', 

/, 5X,  'OF  THE  EXPERIMENTAL  VALUE:  \F6.3, 

'  MILLIMETERS') 

ELSE 

ZDIF2MAX-( ( (Z2+TOLPM)/Z2C)  -  1.)  *  100. 

URITE(* , 212 )Z2CMM , TOLPMMM, Z2HM , ZDIF2MIN , ZDIF2MAX 
FORMAT(// , 5X, 'THE  CALCULATED  VALUE  OF  PIN  T-2:  ' ,F6.3, 

'  MILLIMETERS', 

/.5X, 'IS  NOT  WITHIN  THE  PLACEMENT  TOLERANCE:  \F5.3, 

'  MILLIMETERS'. 

/ ,  5X ,  'OF  THE  EXPERIMENTAL  VALUE:  \F6.3, 

'  MILLIMETERS', 

//,5X,'THE  MINIMUM  PERCENT  DIFFERENCE  BETWEEN  THE  CALCULATED ' , 
/,5X, 'AND  EXPERIMENTAL  VALUES  (CONSIDERING  THE  PLACEMENT', 
/.5X. 'TOLERANCE)  IS:  '.F7.2,'  PERCENT', 

//,5X.'THE  HAXIMUM  PERCENT  DIFFERENCE  BETWEEN  THE  CALCULATED', 
/ , 5X , ' AND  EXPERIMENTAL  VALUES  (CONSIDERING  THE  PLACEMENT', 
/,5X, 'TOLERANCE)  IS:  *,F7.2,'  PERCENT') 

ENDIF 

ELSE  IF( (Z2C- (Z2 -TOLPM) ) ,GT.0)THEN 
WRITE(* , 214 ) Z2CMM , TOLPMMH , Z2HM 
214  FORMAT(//,5X, 'THE  CALCULATED  VALUE  OF  PIN  T-2:  ',F6.3, 

1  '  MILLIMETERS', 

2  /, 5X, 'IS  WITHIN  THE  PLACEMENT  TOLERANCE:  '.F5.3, 

3  '  MILLIMETERS', 

4  /,  5X, 'OF  THE  EXPERIMENTAL  VALUE:  '.F6.3, 

5  '  HILLIMETERS' ) 

ELSEI F(ABS ( Z2C ). LE . ABS ( Z2 • TOLPM) ) THEN 
ZDIF2HIN  -ABS (I  -  ( (Z2-TOLPM)/Z2C) )  *  100 
ZDIF2MAX  -ABS(1  -  ( (Z2+TOLPH)/Z2C) )  *  100 
WRITE (* , 2X6 ) Z2CHM , TOLPMMM , Z2MH , ZDIF2HIN , ZDIF2HAX 
ELSE 

ZDIF2MIN  -ABS( (Z2C/(Z2*TOLPH) )  -  1)  *  100 

IF( ABS ( Z2C ). LE . ABS ( Z2+TOLFH) ) THEN 
ZDIF2MAX  -ABS(1  -  ((Z2+TOLPM)/Z2C))  *  100 
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ELSE 

ZDIF2MAX  -ABS ( (Z2C/(Z2+TOLPM) )  -  1)  *  100 
ENDIF 


WRITE (*,216) Z2CHH , TOLPMMM , Z2MM,  ZDIF2KIN , ZDIF2HAX 
216  FORMAT (// , 5X , ' THE  CALCULATED  VALUE  OF  PIN  T- 2:  \F6.3, 

1  *  MILLIMETERS', 

2  /,5X, 'IS  NOT  WITHIN  THE  PLACEMENT  TOLERANCE:  '.F5.3, 

3  '  MILLIMETERS', 

4  /,5X, ‘OF  THE  EXPERIMENTAL  VALUE:  \F6.3, 

5  1  MILLIMETERS', 

6  //,5X, 'THE  MINIMUM  PERCENT  DIFFERENCE  BETWEEN  THE  CALCULATED', 

7  /, 5X, 'AND  EXPERIMENTAL  VALUES  (CONSIDERING  THE  PLACEMENT*, 

8  /,5X,  'TOLERANCE)  IS:  \F7.2,'  PERCENT', 

9  //,5X, 'THE  MAXIMUM  PERCENT  DIFFERENCE  BETWEEN  THE  CALCULATED', 

1  /, 5X, 'AND  EXPERIMENTAL  VALUES  (CONSIDERING  THE  PLACEMENT', 

2  /,5X, 'TOLERANCE)  IS:  ',F7.2,'  PERCENT') 


ENDIF 


WRITE(*,80) 

READ(*,*) 

C  CHECK  PIN  3  BY  USING  EQN  OF  PLANE  THRU  PINS  1,  2,  4 

COEFF1  -  (Y2-Y1)*(Z4-Z1)  -  <Y4-Y1)*(Z2-ZI) 

COEFF2  -  (X2-X1)*(Z4-Z1)  -  (X4-X1)*(Z2-Z1) 

COEFF3  -  (X2-X1)*(Y4-Y1)  -  (X4-X1)*(Y2-Y1) 

CONST  -  - (X1*C0EFF1)  ♦  (Yl*COEFF2)  -  (Z1*C0EFF3) 
Z3C  -  ((Y3*COEFF2)  -  (X3*C0EFF1)  -  CONST>/COEFF3 
Z3CHH  -  Z3C  *  1000. 


URITE(*. 120)COEFF1 , COEFF2 , COEFF3 . CONST 
120  FORMAT(//,18X, *****  CHECK  LOCATION  OF  TOA  PIN  T-3B  ****', 

1  //, 5X, 'THE  EQUATION  OF  THE  PLANE  THROUGH  TOA  PINS  1,  2,', 

2  '  &  4  IS: ' , 

3  //,10X,EU.S, '  X  +  '  ,E11.5,  *  Y  +  '.Eil.S,'  Z  +  '.E11.5,'  -  O') 
Z3HIN  -  AHIN1(Z3,Z3C) 

IF(Z3MIN.EQ.Z3)  THEN 

C  CALCULATED  VALUE  OF  Z3  IS  GREATER  THAN  THE  EXPERIMENTAL  VALUE  OF  Z3 
ZDIF3HIN  -  ( ( Z3C/ ( Z3+TOLPM ) )  *  1.)  *  100. 

IF(ZDIF3HIN.LE.O. )  THEN 
WRITE (*, 1 22 ) Z3CHM , TOLPMMM , Z3MM 

122  FORMAT(//,5X, 'THE  CALCULATED  VALUE  OF  PIN  T-3B:  '.F6.3, 

1  '  MILLIMETERS ' , 

2  /,$X,'IS  WITHIN  THE  PLACEMENT  TOLERANCE:  \F5.3, 

3  ’  MILLIMETERS', 

4  /,5X, 'OF  THE  EXPERIMENTAL  VALUE:  ',F6.3, 
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'  MILLIMETERS') 


ELSE 

ZDIF3MAX  -  ABS( (Z3C/(Z3-TOLPM) )  -  1.)  *  100. 

WRITE (* , 124) Z3GMM , TOLPMMM , Z3MM, ZDIF3MIN , ZDIF3MAX 
124  FORMAT (//,5X, 'THE  CALCULATED  VALUE  OF  PIN  T-3B:  \F6.3, 

1  '  MILLIMETERS', 

2  / , 5X , ' I S  NOT  WITHIN  THE  PLACEMENT  TOLERANCE :  ' , F5 . 3 , 

3  '  MILLIMETERS', 

4  /,5X, 'OF  THE  EXPERIMENTAL  VALUE:  ',F6.3, 

5  '  MILLIMETERS ' , 

6  //, 5X, 'THE  MINIMUM  PERCENT  DIFFERENCE  BETWEEN  THE  CALCULATED' 

7  /, 5X, 'AND  EXPERIMENTAL  VALUES  (CONSIDERING  THE  PLACEMENT', 

8  /,5X, 'TOLERANCE)  IS:  '.F7.2,'  PERCENT', 

9  // , 5X , 'THE  MAXIMUM  PERCENT  DIFFERENCE  BETWEEN  THE  CALCULATED ' 

1  /,5X, 'AND  EXPERIMENTAL  VALUES  (CONSIDERING  THE  PLACEMENT' , 

2  /,5X, 'TOLERANCE)  IS:  *,F7.2,'  PERCENT') 

ENDIF 

C  CALCULATED  VALUE  OF  Z3  IS  LESS  THAN  THE  EXPERIMENTAL  VALUE  OF  Z3 


ELSE  IF  (Z3C.GT.0)  THEN 
ZDIF3MIN-(((Z3-TOLPH)/Z3C)  -  1.)  *  100. 

IF(ZDIF3MIN.LE.O)  THEN 
WRITE ( * , 2 20 ) Z3CHM , TOLPHHH , Z3MM 

220  FORMAT(//,SX, 'THE  CALCULATED  VALUE  OF  PIN  T-3B:  '.F6.3, 

'  MILLIMETERS ' , 

/,5X,  'IS  WITHIN  THE  PLACEMENT  TOLERANCE:  \F5.3. 

'  MILLIMETERS', 

/ , SX . ' OF  THE  EXPERIMENTAL  VALUE:  \F6.3, 

'  MILLIMETERS') 

ELSE 

ZDIF3MAX-(((23+TOLPM)/Z3C)  -  1.)  *  100. 

WRITE ( * , 2  22 ) Z3CMM , TOLPHHH , Z3MH . ZDIF3MIN , ZDI F3HAX 
222  F0RHAT(//,5X, 'THE  CALCULATED  VALUE  OF  PIN  T-3B:  *,F6.3, 

1  '  MILLIMETERS'. 

2  /,  SX, 'IS  NOT  WITHIN  THE  PLACEMENT  TOLERANCE:  \F5.3, 

3  '  MILLIMETERS', 

4  /.SX.'OF  THE  EXPERI MENTAL  VAUJE :  ' , F6 . 3 , 

5  '  MILLIMETERS'. 

6  //, SX, 'THE  MINIMUM  PERCENT  DIFFERENCE  BETWEEN  THE  CALCULATED' 

7  / , SX , ' AND  EXPERIMENTAL  VALUES  (CONSIDERING  THE  PLACEMENT ' , 

8  /,5X, 'TOLERANCE)  IS:  '.F7.2,'  PERCENT', 

9  //, 5X, 'THE  MAXIMUM  PERCENT  DIFFERENCE  BETWEEN  THE  CALCULATED' 

1  /.sx.’and  experimental  values  (considering  the  placement' , 

2  /,5X, 'TOLERANCE)  IS:  \P7.2,'  PERCENT') 

ENDIF 


ELSE  IF( (Z3C- (Z3-TGLPK) ) .GT. 0)THEN 
WRITE  (* , 224) Z3CHM, TOLPMMM, Z3HH 

224  FORMAT (//, SX, ’THE  CALCULATED  VALUE  OF  PIN  T-3B:  \F6.3, 

1  ’  MILLIMETERS', 

2  /,5X, 'IS  WITHIN  THE  PLACEMENT  TOLERANCE:  '.F5.3. 
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3  ’  MILLIMETERS ' , 

4  /,5X, 'OF  THE  EXPERIMENTAL  VALUE:  \F6.3, 

5  '  MILLIMETERS') 

ELSE  IF(AES(Z3C) . LE.A£S(Z3-TOLPM))THSN 
ZDIF3MIN  ~ABS(1  -  < (Z3-TOLPM)/Z3C))  *  100 
2DIF3MAX  -ABS (1  -  ( (Z3+TOLPM)/Z3C) )  *  100 

WRITE ( * , 226 ) Z3CMH , TOLPMMM , Z3MM , ZDIF3MXN , ZDIF3MAX 
ELSE 

ZDIF3HIN  -ABS ( (Z3C/(Z3 -TOLPM) )  -  1)  *  100 

IF(ABS(Z3C) .LE. ABS (Z3+TOLPM) ) THEN 
ZDIF3MAX  -ABS (I  -  (<Z3+TOLFK)/Z3C))  *  100 
ELSE 

ZDIF3MAX  -AB3( (Z3C/(Z3+TOLPH) )  -  1)  *  100 
ENDIF 


WRITE (*,226 )  23CHM ,  TOLPW&! ,  ?-3MH .  ZDIF3HXN .  ZDIF3SAX 
226  FORMAT (/7.5X, 'THE  CALCULATED  VALUE  OF  MB  T-3B:  \F6.3, 

1  *  MILLIMETERS.  • , 

2  /,5X,  'IS  NOT  WITHIN  THE  PLACEMENT  TOLERANCE:  \P5.3, 

3  '  MILLIMETERS', 

4  / , 5X , ' OF  THE  EXPERIMENTAL  VALUE :  \  F6 . 3 . 

5  '  MILLIMETERS', 

6  //,5X, 'THE  HINIHUM  PERCENT  DIFFERENCE  BETWEEN  THE  CALCULATED’ . 

7  /.SX.'AND  EXPERIMENTAL  VALUES  (CONSIDERING  THE  PLACEMENT1 , 

8  /.SX. 'TOLERANCE)  IS:  \F7.2,'  PERCENT', 

9  //, 5X, 'THE  MAXIMUM  PERCENT  DIFFERENCE  BETWEEN  THE  CALCULATED ' , 

1  /.SX.'AND  EXPERIMENTAL  VALUES  (CONSIDERING  THE  PLACEMENT', 

2  /.SX. 'TOLERANCE)  IS:  \F7.2,*  PERCENT') 


ENDIF 

WRITE (*,80) 

READ(*,*) 

C  CHECK  LOCATION  OF  PIN  4  BY  USING  EQN  OF  PLANE  THRU  PINS  1,  2.  &  3 

CQEFFl  -  (Y2-¥1)*(Z3-Zl)  -  <Y3-Y1)*{Z2-ZI) 

COEFF2  -  (X2-X1)*(Z3-Z1)  *  (X3*X1)*(22*Z1) 

COEFP3  -  (X2*X1)*(Y3-Y1)  -  (X3-XI)*(Y2-Y1) 

CONST  -  - (X1*CQEFF1)  +  (Yl*COEFF2)  -  (ZI*COEFF3) 

Z4C  -  ( (Y4*COEFF2)  -  (X4*CGEFFI)  -  CONST) /COEFF3 
24CHM  -  Z4C  *  1000. 

WRITE(* , 130)C0EFF1 . COEFF2 . COEFF3 , CONST 
130  FORKAT(//,18X, *****  CHECK  LOCATION  OF  TOA  PIN  T-4  ****', 

1  //, SX, 'THE  EQUATION  OF  THE  PLANE  THROUGH  TGA  PINS  1,  2,', 

2  '  &  3  IS:  ' , 

3  //.10X.EU.S, '  X  +  ' , E11.5,  *  Y  +  \E11.5,'  Z  ♦  '.E11.5,'  -  O') 
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Z4MIN  -  AMIN1(Z4,Z4C) 


IF(Z4MIN.EQ.Z4)  THEN 

C  CALCULATED  VALUE  OF  Z4  IS  GREATER  THAN  THE  EXPERIMENTAL  VALUE  OF  Z4 
ZDIF4MIN  -  ((Z4C/(Z4+TOLPM))  -  1.)  *  100. 

IF(ZDIF4MIN.LE.O.)  THEN 
WRITE ( * , 1 3  2 ) Z4CMM , TOLPMMM , Z4MM 

132  FORMAT(//,5X, 'THE  CALCULATED  VALUE  OF  PIN  T-4:  \P6.3, 

1  '  MILLIMETERS', 

2  /, 5X, ' IS  WITHIN  THE  PLACEMENT  TOLERANCE:  '.F5.3, 

3  '  MILLIMETERS ' , 

4  / , 3X , s OF  THE  EXPERIMENTAL  VALUE :  ' , F6 . 3 , 

5  '  MILLIMETERS') 

ELSE 

ZDIF4MAX  -  ((Z4C./(Z4-TOLPM) )  -  1.)  *  100. 


WRITE {* , 134 ) Z4CMM . TOLPHMM . 24MM , ZDIF4HIN , ZUIF4HAX 
134  FORMAT <// ,  5X , ' THE  CALCULATED  VALUE  OF  PIN  T-4:  *,F6.3, 

1  '  MILLIMETERS', 

2  / , 5X, 1  IS  NOT  WITHIN  THE  PLACEMENT  TOLERANCE:  \F5.3, 

3  '  MILLIMETERS', 

4  / , 5X . ! OF  THE  EXPERIMENTAL  VALUE :  ' , F6 . 3 , 

5  '  MILLIMETERS', 

f-  // ,  5X ,  *  THE  MINIMUM  PERCENT  DIFFERENCE  BETWEEN  THE  CALCULATED', 

7  /,5X, 'AND  EXPERIMENTAL  VALUES  (CONSIDERING  THE  PLACEMENT', 

8  /,5X, 'TOLERANCE)  IS:  \F7.2,'  ''ERCENT ’  , 

9  //, 5X, 'THE  MAXIMUM  PERCENT  DIFFERENCE  BETWEEN  THE  CALCULATED', 

1  /, 5X, 'AND  EXPERIMENTAL  VALUES  (CONSIDERING  THE  PLACEMENT* , 

2  /,5X, 'TOLERANCE)  IS:  \F7.2,*  PERCENT') 

ENDIF 


C  CALCULATED  VALUE  OF  Z4  IS  LESS  THAN  THE  EXPERIMENTAL  VALUE  OF  Z4 


ELSE  IF  (Z4C.CT.0)  THEN 
ZDIP4HIN-(((Z4-T0LPH)/Z4C)  -  1.)  *  100. 

IF(2DIF4HIN.l£.G)  THEN 
WRITE (* ,230) 24CHM . TOLPMMM . Z4HH 

230  FORMAT (//,  SX. 'THE  CALCULATED  VALUE  OF  PIN  T'4:  \F6.3. 

1  '  MILLIMETERS', 

2  /.SX.'IS  WITHIN  THE  PLACEMENT  TOLERANCE:  \F5.3, 

3  '  MILLIMETERS', 

4  / , 5X , ' OF  THE  EXPER) MENTAL  VALUE :  ' , F6 . 3 , 

5  '  MILLIMETERS') 

ELSE 

ZDIF4MAX  -  ABS( ( (Z4+TOLPH)/Z4C)  -  1.)  *  100. 

WR ITE  (* ,  2  32 )  Z4CHM ,  TOLPMMM ,  Z4MM ,  ZD!  F4MIN .  ZDIF4HAX 
232  FORMAT(//, 5X, ‘THE  CALCULATED  VALUE  OF  PIN  T-1B:  ',F6.3, 

1  *  MILLIMETERS', 

2  /.SX.'IS  NOT  WITHIN  THE  PLACEMENT  TOLERANCE:  \F5.3, 

3  '  NILLIHETERS ' , 

4  / , SX , *  OF  THE  EXPERIMENTAL  VALUE :  * . F6 . 3 , 

3  '  MILLIMETERS', 

6  //.SX.'THE  HINIHUH  PERCENT  DIFFERENCE  BETWEEN  THE  CALCULATED*, 


E- 18 


7  /, 5X, 'AND  EXPERIMENTAL  VALUES  (CONSIDERING  THE  PLACEMENT', 

8  /,5X, 'TOLERANCE)  IS:  '.F7.2,'  PERCENT', 

9  //,5X, 'THE  MAXIMUM  PERCENT  DIFFERENCE  BETWEEN  THE  CALCULATED', 

1  /,5X, 'AND  EXPERIMENTAL  VALUES  (CONSIDERING  THE  PLACEMENT' , 

2  /,5X, 'TOLERANCE)  IS:  ',F7.2,'  PERCENT') 

ENDIF 

ELSE  IF( (Z4C- (Z4-T0LPM) ) .GT.0)THEN 
WRITE (*, 234) Z4CMM , TOLPMMM , Z4MM 
FORMAT (  // ,  5X",  '  THE  CALCULATED  VALUE  OF  PIN  T-4:  '.F6.3, 

'  MILLIMETERS’, 

/.SX.'IS  WITHIN  THE  PLACEMENT  TOLERANCE:  \F5.3, 

'  MILLIMETERS', 

/,5X, 'OF  THE  EXPERIMENTAL  VALUE:  *,F6.3, 

'  MILLIMETERS') 

ELSE  XF(ABS(Z4C).LE.ABS(Z4-TOLPM))THEN 
ZDIF4MIN  “AB3  ( 1  -  ( (Z4-TOLPMVZ4C) )  •*  100 
ZDIF4MAX  ~ABS(1  -  ( (Z4+TOLPM)/Z4C) )  *  100 

WRITE(* , 2  36) Z4CMM , TOLPMMM , Z4MM , ZDIF4MIN , ZDIF4MAX 
ELSE 

ZDIF4MIN  -ABS( (Z4C/(Z4-T0LPM) )  -  I)  *  100 

IF(ABS(?4C) . LE . ABS ( Z4rTOLPM) ) THEN 
ZDIF4NAX  -ABS(1  *  ((Z4+TOLPM)/Z4C))  *  100 
ELSE 

ZDIF4MAX  -ABS ( ( Z4C/ ( Z4+TOLPM) )  -  1)  *  100 
ENBIF 
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WRITE ( * , 2  36 ) Z4CMK , TOLPMMM , Z4MM , ZDI F4HIN , ZDIF4MAX 
236  FORMAT (//, 5X, ’THE  CALCULATED  VALUE  OF  PIN  T-4:  '.F6.3, 

1  '  MILLIMETERS', 

2  /,5X,' IS  NOT  WITHIN  THE  PLACEMENT  TOLERANCE:  V.F5.3, 

3  '  MILLIMETERS ’, 

4  / , 5X , ' OF  THE  EXPERIMENTAL  VALUE :  ' , F6 . 3 , 

5  1  MILLIMETERS', 

6  //, 5X, 'THE  MINIMUM  PERCENT  DIFFERENCE  BETWEEN  THE  CALCULATED' , 

7  /. 5X, ’AND  EXPERIMENTAL  VALUES  (CONSIDERING  THE  PLACEMENT', 

8  /,5X, 'TOLERANCE)  IS:  '.F7.2,'  PERCENT' , 

9  // , SX , ' THE  MAXIMUM  PERCENT  DIFFERENCE  BETWEEN  THE  CALCULATED', 

1  /,5X,  'AND  EXPERIMENTAL  VALUES  (CONSIDERING  THE  PLACEMENT', 

2  /,5X, 'TOLERANCE)  IS:  '.F7.2,'  PERCENT') 


ENDIF 


WRITE(*,80) 

R£AD(*,*> 

1000  END 
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